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SUMMARY 


The  standard  procedure  used  in  estimating  the  source  size  of  underground  nuclear 
explosions  using  mi,  measurements  has  been  to  separate  the  station  terms  from  the 
network-averaged  source  terms.  The  station  terms  thus  derived  actually  reflect  the  combina¬ 
tion  of  the  path  effect  and  the  station  effect,  when  only  those  events  in  a  close  proximity  are 
utilized.  If  worldwide  explosions  are  used  in  the  inversion,  then  the  path  effect  tends  to  be 
averaged  out  at  each  station.  In  either  case,  the  effect  due  to  the  propagation  path  alone 
would  not  be  obvious. 

Under  this  research  contract  (F29601-91-C-DB23),  we  propose  to  decompose  the  station 
amplification  effect  further  with  a  joint  inversion  scheme  which  simultaneously  determines  the 
seismic  source  size,  the  path  terms,  and  the  receiver  terms.  Short-period  P-wave  amplitudes 
of  217  worldwide  underground  nuclear  explosions,  including  92  blasts  from  Semipalatinsk, 
recorded  at  118  WWSSN  stations  have  been  used  in  one  single  inversion  to  infer  the  2733 
unknown  parameters.  For  Semipalatinsk  events  in  our  WWSSN  database,  the  new  factor¬ 
ing  procedure  provides  more  stable  measurements  across  the  whole  recording  network 
with  a  reduction  in  the  fluctuational  variation  by  a  factor  of  up  to  3.  A  list  of  station  corrections 
(which  are  applicable  to  all  nuclear  test  sites  around  the  world)  and  the  path  corrections  for 
five  Central  Asian  test  sites  is  compiled.  If  WWSSN  recordings  of  a  new  event  become  avail¬ 
able,  these  corrections  can  be  applied  immediately  without  re-running  the  inversion  in  a  batch 
mode. 

We  have  recomputed  the  yield  estimates  of  92  Semipalatinsk  explosions  based  on  the 
path-corrected  values  derived  in  this  study.  Based  on  the  first  motion  of  the  P  wave 
alone,  the  (central-value)  yield  estimate  of  Soviet  JVE  (09/14/88)  is  113  KT  and  that  of  the 
cratering  event  on  01/15/65  at  Balapan  test  site  is  111  KT,  which  are  in  excellent  agreement 
with  estimates  derived  by  other  means.  Thus  the  first  motion  of  the  initial  short-period  P 
waves  appears  to  be  a  very  favorable  source  measure  for  explosions  fired  in  hard  rock  sites 
underlain  by  the  stable  mantle  (such  as  Semipalatinsk).  The  bias  relative  to  NTS  at  50KT 
level  is  inferred  as  0.36  magnitude  unit. 

Also  included  in  this  Semi-annual  Report  No.  /  is  a  preliminary  assessment  of  WWSSN’s 
capability  in  remotely  monitoring  Eurasian  explosions.  A  strong  correlation  between  the 
P-wave  amplitude  and  Lg  detection  at  teleseismic  distance  is  observed.  Further  study  along 
this  line  is  needed  to  investigate  its  implication  to  low-yield  threshold  global  monitoring. 
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1.  tNTRODUCTION 


The  main  problem  with  the  conventional  is  that  it  is  a  rather  nebulous  parameter; 
simply,  it  is  a  function  of  the  largest  peak-to-peak  amplitude  in  the  first  few  seconds  of  P 
wave  motion  with  adjustment  for  the  period  of  the  arriving  phase.  The  parameter  mt,  was 
adapted  from  the  need  to  order  systematically  the  size  of  earthquakes.  The  measure  itself  has 
inherent  impreciseness  as  the  measure  is  not  related  to  the  physics  of  the  source  per  se  but 
is  the  largest  constructive  interference  of  waves  originating  at  the  source,  source  region,  pro¬ 
pagation  path,  receiver  region,  and  receivers  (Butler,  1981;  Johnson,  1981).  To  relate  the 
to  the  seismic  yield,  all  effects  not  due  to  the  source  must  naturally  be  corrected.  It  is  often 
difficult,  however,  to  separate  these  effects.  In  fact,  the  effects  of  source  and  propagation  are 
often  indistinguishable,  and  unless  one  is  known  the  other  cannot  be  uniquely  determined 
(Johnson,  1981).  Consequently,  it  was  reported  to  be  difficult  to  make  measurements  that 
are  internally  consistent  within  0.1  mt,  with  the  conventional  mi,  (Bache,  1982),  simply 
because  some  of  the  aforementioned  effects  are  not  accounted  for  accurately. 

Von  Seggem  (1973)  showed  that  including  station  corrections  typically  halves  the  stan¬ 
dard  deviation  of  mt,  from  North  American  LRSM  stations  recording  NTS  events.  Even  better 
results,  with  the  standard  deviation  reduced  by  a  factor  of  3  or  so,  can  be  obtained  for  a  net¬ 
work  with  all  stations  beyond  30®.  Applying  station  corrections  to  the  determination  or  the 
network  spectra  averaging  has  become  a  standard  procedure  in  this  nuclear  monitoring  com¬ 
munity  {e.g..  Ulwall  etal.,  1988;  Murphy  etal..  1989;  Sykes  and  Ekstrom,  1989;  Jih  and 
Shumway,  1989;  and  many  others).  The  station  effects  are  strongly  dependent  on  azimuth 
(Chang  and  von  Seggern,  1980),  which  led  Bache  (1982)  and  many  others  to  believe  that  sta¬ 
tistical  station  corrections  will  not  be  nearly  so  effective  in  reducing  m^  variance  in  the  multi¬ 
ple  source  region  problem. 

Marshall  etal.  (1979)  attempted  to  correct  several  important  factors  that  can  bias  mi,. 
They  used  bulletin  log(A/T)  data.  These  data  are  corrected  for  receiver-station  attenuation 
differences,  and  the  resulting  magnitude  is  called  m^.  Correcting  m2  for  source-region 
attenuation  gives  m3,  and  correcting  m3  for  source  depth  gives  mQ.  (The  network  averages 
of  these  m^s  are  denoted  by  ,  mj,  m3,  and  mQ,  respectively.)  The  major  change  in  this 
scheme  is  associated  with  the  source-region  correction,  which  can  be  as  different  as  0.4  m.u. 
between  sites  or  a  factor  of  about  2.5  in  yield  estimates.  Essentially  this  approach  is  based 
on  the  discovery  by  Marshall  and  Springer  (1976)  and  Douglas  etal.  (1981)  that  LRSM 
amplitude  resiauals  correlate  with  P„  velocity  near  the  stations  and  the  assumption  that  such 
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correlation  between  the  attenuation  and  Is  valid  elsewhere  as  well.  However,  it  turns  out 
thd*  the  standard  deviation  of  the  mo  is  not  less  than  that  for  the  m,  and  m2  from  the  same 
data  set  (Marshall  etai,  1979;  Bache,  1982). 

It  is  obvious  that  the  only  way  to  reduce  the  statistical  fluctuation  is  to  obtain  fundamental 
causal  knowledge  of  the  focusing  and  defocusing  beneath  the  source  and  receiver.  We 
expect  teleseismic  P-wave  amplitudes  to  vary  as  the  source  location  changes  within  a  test 
site.  The  mu  residuals  (with  respect  to  the  bast-fitting  -yield  curve)  of  NTS  events  show 
systematic  trends  that  are  consistent  with  local  tectonics  (Minster  et  at. ,  1981).  At  Yucca  Flat, 
the  residuals  are  positive  to  the  west  and  negative  to  the  east  of  the  north-south  trending  nor¬ 
mal  fault  system  that  bisects  the  valley.  At  Pahute  Mesa,  the  spatial  pattern  is  less  clear,  but 
the  residuals  tend  to  be  negative  toward  the  center  of  the  buried  Silent  Canyon  Caldera  and 
positive  toward  the  edges.  An  attractive  explanation  is  that  these  variations  are  due  to 
focusing/defocusing  effects  that  are  not  averaged  out  over  the  network,  although  the  possibil¬ 
ity  of  systematic  source-coupling  difference  has  not  been  eliminated. 

Jih  and  Wagner  (1991b,  1992)  propose  to  compute  the  new  station  magnitude  012  9  for 
the  i-th  event  recorded  at  the  j-th  station  as 

^2,9(1.])  =  logio[A(i,j)/T{i.j)]  +  B(A(i,j))  -  S(i)  -  F(k(i),j) 

where  A(i,j)  is  the  displacement  amplitude  (in  millimicrons)  and  T(i,j)  is  the  period  (in  seconds) 
of  the  P  wave.  The  B(A)  is  the  distance-correction  term.  S(j)  is  the  station  correction,  and 
F(k(i),j)  is  the  path  correction  for  explosions  from  the  k(i)-th  source  region.  The  resulting  new 
magnitude  is  called  012  9  to  avoid  confusion  with  ttie  m3  defined  in  Marshall  et  at.  (1979)  that 
corrects  for  the  source-region  attenuation  and  station  terms  solely  based  on  published  P„ 
velocity.  The  path  corrections  determined  in  this  procedure  provide  direct  and  informative 
clues  to  characterize  the  various  propagation  paths.  We  will  also  use  m29  extensively 
throughout  this  study  to  characterize  the  magnitude-yield  scaling  relationship  for  Semipala- 
tinsk  explosions. 
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2.  PATH-CORRECTED  UNBIASED  NETWORK  ESTIMATOR 

The  conventional  definition  of  the  station  magnitude  is  computed  as 

^6  =  logio(A/T)  +  B(A)  [1] 

where  A  is  the  displacement  amplitude  (in  nm)  and  T  is  the  predominant  period  (in  sec)  of  the 
P  wave.  The  B(A)  is  the  distance-correction  term  that  compensates  for  the  change  of  P- 
wave  amplitudes  with  distance  (e.g..  Gutenberg  and  Richter,  1956;  Veith  and  Clawson,  1972). 
nrii,  in  [1]  is  also  denoted  as  m,  in  Marshall  et  al.  (1979).  The  ISC  bulletin  rrib  is  just  the  net¬ 
work  average  of  these  raw  station  nib  values  without  any  further  adjustment. 

Consider  Ng  explosions  detonated  at  Np  source  regions  that  are  recorded  at  some  or  all 
of  Ns  stations.  The  conventional  GLM  [General  Linear  Model]  or  LSMF  [Least  Squares  Matrix 
Factorization]  network  (Douglas,  1966;  Blandford  and  Shumway,  1982;  Marshall 
efa/.  ,1984;  Jih  and  Shumway,  1989;  Murphy  etai,  1989)  is  the  least-squares  or  maximum- 
likelihood  network  average  of  the  “station-corrected"  magnitudes; 

/772.2(M)  =  "Ji(M>-S0)  [2] 

where  SG)  is  the  “statistical"  or  “empirical"  receiver  correction  at  the  j-th  station.  In  Marshall 
efa/.  (1979),  a  priori  information  about  the  P„  velocity  underneath  each  station  is  used  to 
determine  its  associated  “deterministic"  receiver  correction,  S(j),  and  the  resulting  magnitude 
is  called  m2.  The  GLM  receiver  corrections,  however,  are  inferred  jointly  from  a  suite  of 
event-station  pairs,  and  no  a  priori  geophysical  or  geological  condition  is  assumed  (and 
hence  the  different  notation  m2,2).  The  high  correlation  between  the  tectonic  type  and  the 
GLM  station  terms  suggests  that  the  empirical  station  corrections  do  reflect  the  averaged 
upper  mantle  conditions  underneath  the  receivers,  if  the  azimuthal  coverage  at  each  station  is 
broad  enough. 

Jih  and  Wagner  (1991b,  1992)  propose  to  compute  the  new  station  magnitude  ^39  for 
the  i-th  event  recorded  r  the  j-th  station  as 

rn2,9(i,j)  =  m,  (i,j)  -  SG)  -  F(k(i),j)  =  mj  2  (i.j)  -  F(k(i),j)  [3] 

At  the  j-th  station,  F(k(*),j)  is  a  constant  for  all  events  detonated  in  the  same  k-th  “geologically 
and  geophysically  uniform  region".  Partitioning  a  sino'e  nuclear  test  site  into  several 
“regions”  may  be  necessary  in  order  to  account  accurately  for  the  focusing/defocusing  effects. 
This  012  9  is  very  similar  to  the  m^  in  Marshall  et  al.  (1979)  except  that,  again,  a  priori 
attenuation  information  of  the  source  region  is  used  in  Marshall  et  al.  (1979)  to  determine  the 


3 


correction  term,  whereas  Jih  and  Wagner  (1991a,  1991b.  1992)  invert  for  the  path  or  near¬ 
source  effects  from  the  data  empirically.  In  other  words,  the  source-region  corrections  pro¬ 
posed  by  Marshall  eta/.  (1979)  are  constants  (for  all  explosions  in  the  same  source  region) 
regardless  of  the  location  of  the  seismic  stations,  whereas  the  path  or  near-source  corrections 
in  Equation  [3]  are  highly  dependent  on  the  source-station  paths. 

We  now  examine  briefly  the  fundamental  difference  between  the  present  scheme  (Equa¬ 
tion  [3])  and  the  previous  GLM  schemes.  In  LSMF  and  the  standard  GLM  schemes  (e.g., 
Douglas,  1966;  Blandford  and  Shumway,  1982;  Marshall  etai,  1984;  Lilwall  eiai,  1988;  Jih 
and  Shumway,  1989;  Murphy  etai,  1989),  it  is  assumed  that  the  observed  station  mM  is 
the  sum  of  the  true  source  size  of  the  i-th  event,  E(i),  the  receiver  term  of  the  j-th  station,  S0, 
and  the  random  noise,  v(i,j): 

mi,(i.j)  =  E(i)-hSG)  +  v(i,j)  [4] 

The  receiver  term,  SG).  is  constant  with  respect  to  all  explosions  from  different  test  sites,  and 
hence  it  would  inherently  reflect  the  “averaged”  receiver  effect  —  provided  the  paths  reaching 
the  station  have  broad  azimuthal  coverage.  When  world-wide  explosions  are  used,  the  stan¬ 
dard  deviation  (a)  of  the  noise  v  in  [4]  is  typically  about  0.3  m.u.  or  larger. 

If  LSMF  or  GLM  is  applied  to  events  within  a  smaller  area  of  source  region,  then  the  a  of 
V  in  [4]  could  reduce  to  0.15-0.2  m.u.  However,  the  result  of  such  “single-test-site  GLM” 
approach  should  be  interpreted  or  utilized  cautiously.  The  event  m,,  values  (/.e. ,  the  "E"  term 
in  [4])  so  determined  are  excellent  estimates  of  “relative  source  size"  for  that  test  site  only.  If 
this  “single-test-site  GLM"  inversion  is  applied  to  several  test  sites  separately,  it  may  not  be 
easy  or  obvious  to  find  a  consistent  baseline  for  “absolute  yield”  estimation  or  immediate 
combination  of  the  (inter-site)  magnitudes,  since  the  recording  network  is  typically  different 
from  one  test  site  to  another,  and  hence  the  station  terms  are  inevitably  inconsistent.  Further¬ 
more,  the  station  terms  derived  by  the  “single-test-site  GLM”  may  not  necessarily  represent 
the  attenuation  underneath  the  receiver  side  alone.  They  could  be  “contaminated"  or  some¬ 
times  even  overwhelmed  by  the  path/near-source  effects  shared  by  the  explosions  confined  in 
a  narrow  azimuthal  range.  This  could  explain  the  once  puzzling  and  controversial 
phenomenon  Butler  and  Ruff  (1980)  (also  Butler,  1981;  Burdick,  1981)  reported,  namely  that 
using  Soviet  explosions  from  one  test  site  alone  may  fail  to  discern  the  attenuation  differential 
between  the  eastern  and  western  U.S.  There  is  no  doubt,  however,  that  the  GLM  or  LSMF 
type  of  methodology  can  infer  the  station  terms  which  are  strongly  correlated  with  the  upper 
mantle  attenuation  underneath  the  stations,  provided  the  seismic  sources  have  a  broad  spatial 
coverage  as  did  those  in  North  (1977),  Douglas  and  Marshall  (1983),  Lilwall  and  Neary 
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(1985),  Ringdal  (1986),  Jih  and  Wagner  (1991b.  1992),  and  many  others. 

In  the  present  scheme  ([3]),  however,  we  reformulate  the  whole  model  as 

"Jb(i.j)  =  E(i)  +  ?'j)  +  F{k(i).j)  +  v(i,j)  [5] 

where  F(k(i),j)  is  the  correction  term  at  the  j-th  station  for  the  propagation  effect  or  the  near¬ 
source  focusing/defocusing  effect,  which  is  constant  for  all  events  in  the  k-th  “geologicaliy  and 
geophysically  uniform  region".  For  each  seismic  station,  this  F  can  be  regarded  as  its  azimu¬ 
thal  variation  around  the  mean  station  term  S.  However,  as  explained  previously,  it  would  be 
more  appropriate  to  consider  F  the  path  or  near-source  term  because  the  back  azimuths  at 
the  station  could  be  nearly  identical  for  adjacent  test  sites  (such  as  Degelen  and  Murzhik), 
and  yet  the  "F"  terms  could  be  very  different.  By  incorporating  the  F  term  into  the  model,  the 
a  for  world-wide  explosions  is  reduced  to  about  0.2,  roughly  the  same  level  that  which  a 
"single-test-site  GLM”  could  achieve.  Intuitively,  the  present  scheme  (Equations  [3]  or  [5]) 
provides  a  more  detailed  (and  hence  better)  model  than  that  of  Equation  [4]  in  describing  the 
whole  propagation  path  from  the  source  towards  the  receiver.  Simply  put.  Equation  [4]  yields 
a  stronger  fluctuation  in  the  source  terms,  E,  as  well  as  a  larger  standard  deviation  of  v 
because  each  term  in  the  right-hand  side  of  Equation  [4]  would  have  to  “absorb”  part  of  the 
missing  F  term  in  [5].  This  is  exactly  the  same  reason  why  1x12  2  has  smaller  variation  than 
m,. 

Roughly  speaking,  the  procedure  described  in  [3]  has  the  following  advantages: 

•  It  provides  more  stable  mb  measurements  across  the  whole  recording  network,  as  com¬ 
pared  to  the  conventional  GLM  or  LSMF  procedure  which  only  corrects  for  the  station 
terms.  The  reduction  in  the  standard  deviation  of  network  from  to  ^29  could 
reach  a  factor  of  nearly  3.  As  a  result,  the  scatter  in  ^29  versus  log(yield)  is  smaller 
than  that  for  other  mb- 

•  The  resulting  network  mb  values  are  not  significantly  different  from  the  GLM  results. 
Thus  if  the  network  mb  values  derived  by  GLM  or  LSMF  are  unbiased,  so  are  the  refined 
results. 

•  The  separation  of  the  path  effect  from  the  station  effect  is  a  crucial  step  to  study  the  vari¬ 
ous  propagation  phenomena,  which  in  turn  would  improve  our  understanding  of  the 
seismic  source  as  well. 

We  have  applied  this  procedure  to  217  worldwide  explosions,  and  the  resulting  m2  9 
values  of  92  Semipalatinsk  explosions  are  listed  in  Table  1.  The  118  WWSSN  [World  Wide 
Standarized  Seismograph  Network]  stations  are  selected  such  that  each  station  records  1 0  or 
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more  good  explosion  signals.  There  are  13840  signals,  9080  noise  measurements,  and  1609 
clips  from  17  test  sites  that  are  used  to  invert  for  the  2733  unknown  parameters  with  the 
maximum-likelihood  approach.  The  standard  deviation  of  v{i,j)  in  [5]  is  0.196,  as  compared  to 
0.294  if  the  conventionai  GLM  (Equation  [4])  is  applied  to  the  same  data  set.  A  list  of  the 
event  magnitudes  and  the  path  corrections  for  northern  test  site  in  Novaya  Zemiya  has  been 
presented  in  Jih  and  Wagner  (1992).  Here  we  limit  the  discussion  to  Semipalatinsk  explosions 
only.  The  mi,{Pmax)  25  Balapan  events  for  which  the  mb(Pb)  and  /T7/,(Pa)  are  missing  are 
based  on  station  values  published  by  Uiwall  ef  a/.  (1988).  For  these  station  mt,  values,  a 
compensating  correction  to  convert  the  B(A)  {cf.  Equation  [1])  of  Gutenberg  and  Richter 
(1956)  to  that  of  Veith  and  Clawson  (1972)  is  applied  to  every  station  before  these  record¬ 
ings  are  incorporated  into  our  data  set. 
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Table  1.  Magnitudes  of  Semipalatinsk  Explosions 


Event 

#  of  Signals 

Magnitudes  [7^2.9] 

Yield 

Date 

Site’ 

Ns  Nn  Nc^ 

S.E.M.^ 

P. 

Pb 

p 

'^max 

6501 15B 

BTZ 

45  1  2 

0.028 

5.473 

5.709 

5.865 

100-150 

651121D 

Deg 

48  15  1 

0.024 

4.962 

5.240 

5.452 

29 

66021 3D 

Deg 

51  4  10 

0.024 

5.717 

5.965 

6.152 

125 

660320D 

Deg 

49  9  8 

0.024 

5.416 

5.697 

5.916 

100 

660507D 

Deg 

9  26  1 

0.033 

4.089 

4.237 

4.529 

4 

661019D 

Deg 

51  10  5 

0.024 

5.164 

5.423 

5.596 

20-150 

661218M 

Mzk 

55  8  1 

0.024 

5.395 

5.632 

5.852 

20-150 

670226D 

Deg 

48  9  6 

0.025 

5.438 

5.688 

5.914 

20-150 

6709 16M 

Mzk 

36  29  2 

0.024 

4.657 

4.937 

5.182 

<20 

670922M 

Mzk 

35  31  1 

0.024 

4.516 

4.840 

5.118 

10 

671122M 

Mzk 

7  63  0 

0.023 

3.975 

4.353 

<20 

68061 9B 

BNE 

28  3  2 

0.034 

4.666 

5.002 

5.256 

<20 

680929D 

Deg 

50  7  6 

0.025 

5.222 

5.511 

5.710 

60 

690531 M 

Mzk 

30  30  0 

0.025 

4.468 

4.885 

5.115 

<20 

690723D 

Deg 

38  20  1 

0.025 

4.711 

5.022 

5.248 

16 

69091  ID 

Deg 

19  38  0 

0.026 

4.141 

4.381 

4.709 

<20 

691130B 

BTZ 

49  0  0 

0.028 

5.362 

5.733 

5.915 

125 

691228M 

Mzk 

45  9  3 

0.026 

5.264 

5.551 

5.753 

46 

700721 M 

Mzk 

38  20  1 

0.025 

4.689 

5.033 

5.281 

<20 

701104M 

Mzk 

38  22  1 

0.025 

4.934 

5.137 

5.349 

<20 

710322D 

Deg 

43  14  3 

0.025 

5.117 

5.408 

5.587 

20-150 

710425D 

Deg 

37  5  0 

0.030 

5.434 

5.696 

5.891 

90 

710606M 

Mzk 

38  12  2 

0.027 

4.879 

5.218 

_ 

5.425 

16 

1)  BSW  =  SW  subsite,  Balapan;  BNE  =  NE  subsite,  Balapan;  BTZ  =  transition  zone.  Balapan;  Deg  =  Degelen  Moun¬ 
tain;  Mzk  =  Murzhik. 

2)  Ns  =  #  of  signals,  Nn  =  #  of  noise  measurements.  No  =  <f  of  clips. 


3)  standard  error  in  the  mean. 
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Table  1.  Magnitudes  of  Semipalatinsk  Explosions  (continued) 


Magnitudes  [^2,9] 


Event 

#  of  Signals 

Date 

Site 

Ns  Nn  Nc 

S.E.M. 

710619M 

Mzk 

41  13  0 

0.027 

710630B 

BTZ 

31  19  1 

0.027 

711009M 

Mzk 

27  12  3 

0.030 

711 021 M 

Mzk 

32  9  0 

0.031 

711230D 

Deg 

16  3  0 

0.045 

72021 OB 

BNE 

34  8  2 

0.029 

720328D 

Deg 

28  17  0 

0.029 

72081 6D 

Deg 

23  23  1 

0.029 

720826M 

Mzk 

29  15  2 

0.029 

720902M 

Mzk 

15  29  0 

0.029 

721102B 

BSW 

42  1  15 

0.026 

721210D 

Deg 

30  7  5 

0.030 

721 21 OB 

BNE 

44  2  11 

0.026 

730723B 

BTZ 

52  1  1 

0.027 

731214B 

BNE 

49  8  6 

0.025 

750427B 

BNE 

18  1  1 

0.044 

760704B 

BTZ 

38  0  5 

0.030 

761 123B 

BNE 

22  0  0 

0.042 

761207B 

BSW 

17  21 

0.044 

770329D 

Deg 

25  14  0 

0.031 

770730D 

Deg 

21  16  0 

0.032 

780326D 

Deg 

25  6  0 

0.035 

780422D 

Deg 

21  9  0 

0.036 

4.863 


4.448 


4.791 


4.875 


5.080 


4.811 


4.481 


4.447 


4.688 


4.148 


5.619 


5.075 


5.743 


5.248 


4.904 


5.199 


4.928 


4.401 


4.296 


4.995 


4.562 


5.162 


4.766 


5.026 


5.208 


5.425 


5.073 


4.826 


4.735 


5.033 


4.405 


5.935 


5.402 


5.801 


5.985 


5.549 


5.238 


5.545 


5.351 


4.785 


4.692 


5.301 


4.821 


5.392 


5.036 


5.226 


5.442 


5.610 


5.306 


5.051 


4.991 


5.258 


4.682 


6.158 


5.624 


5.998 


6.171 


5.760 


5.521 


5.812 


5.680 


5.581 


5.073 


4.943 


5.530 


5.071 


2 


16 


6 


8 


<20 


2 


165 


20-150 


140 


Event 


Table  1.  Magnitudes  of  Semipalatinsk 


#  of  Signals 


Explosions  (continued) 


Magnitudes  [^2.9] 


78061 IB 


780705B 


780728D 


780829B 


780915B 


781104B 


781129B 


790623B 


790804B 


790818B 


791028B 


791202B 


791223B 


800522D 


800629B 


800914B 


801012B 


801214B 


801227B 


810422B 


810913B 


811018B 


Site 

Ns  Nn  No 

S.E.M. 

BSW 

170  1 

0.046 

BSW 

38  7  7 

0.027 

Deg 

36  9  6 

0.027 

BNE 

1600 

0.049 

BT2 

36  1  6 

0.030 

BNE 

40  9  6 

0.026 

BSW 

28  0  0 

0.037 

BT2 

40  2  3 

0.029 

BNE 

30  0  0 

0.036 

BSW 

40  4  20 

0.024 

BNE 

28  0  0 

0.037 

BNE 

445  13 

0.025 

BSW 

1500 

0.050 

BSW 

40  3  17 

0.025 

Deg 

36  22  1 

0.025 

BSW 

46  5  6 

0.026 

BT2 

34  5  6 

0.029 

BNE 

23  0  0 

0.041 

BT2 

28  0  0 

0.037 

BNE 

24  0  0 

0.040 

BSW 

25  0  0 

0.039 

BT2 

1700 

0.047 

BSW 

41  3  7 

0.027 

5.246 


5.215 


5.068 


5.414 


5.109 


5.639 


5.609 


5.463 


5.599 


4.721 


5.202 


5.493 


5.513 


5.489 


5.365 


5.655 


5.349 


5.878 


5.894 


5.700 


5.890 


4.980 


5.455 


5.824 


5.492  5.778 


5.811 


5.738 


5.577 


5.926 


5.828 


5.566 


5.886 


6.084 


5.812 


6.114 


6.095 


5.932 


5.900 


6.139 


5.188 


5.664 


6.087 


5.856 


5.919 


5.899 


5.922 


6.077 


5.989 
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Table  1.  Magnitudes  of  Semipalatinsk  Explosions  (continued) 


Magnitudes 


Event 

#  of  Signals 

Date 

Site 

Ns  Nn  Nc 

S.E.M. 

811129B 

BSW 

37  11  5 

0.027 

811227B 

BSW 

23  0  0 

0.041 

820425B 

BTZ 

14  0  0 

0.052 

820704B 

BTZ 

21  0  0 

0.043 

820831 B 

BSW 

27  17  1 

0.029 

821205B 

BSW 

26  0  0 

0.038 

821226B 

BNE 

38  10  1 

0.028 

8306 12B 

BTZ 

1600 

0.049 

831006B 

BSW 

25  0  0 

0.039 

831026B 

BTZ 

18  0  0 

0.046 

831120B 

BNE 

17  83 

0.037 

840425B 

BTZ 

21  0  0 

0.043 

840526B 

BNE 

31  0  3 

0.034 

840714B 

BTZ 

22  0  0 

0.042 

841027B 

BSW 

1900 

0.045 

841202B 

BNE 

22  0  0 

0.042 

841216B 

BTZ 

15  0  0 

0.050 

841228B 

BSW 

19  0  0 

0.045 

85021 OB 

BSW 

18  1  4 

0.041 

850615B 

BSW 

15  0  0 

0.050 

850630B 

BSW 

37  3  6 

0.029 

870620B 

BSW 

24  3  13 

0.031 

8809 14B 

BSW 

25  0  1 

0.038 

5.044  5.313 


4.559  4.865 


5.171  5.378 


4.933  5.130 


5.547  5.848 


5.406 
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3.  RECEIVER  AND  PATH  EFFECTS  ON  P  WAVES  FROM  SEMIPALATINSK 


Figure  1  shows  our  receiver  terms  which  are  inferred  jointly  along  with  the  source-size 
estimates  and  path  terms  from  the  worldwide  explosions  (top)  along  with  the  the  deterministic 
station  terms  predicted  by  Marshall  etal.  (1979)  based  on  published  P„  velocity  around  the 
world  (bottom).  The  receiver  corrections  derived  with  our  approach  match  the  average  tec¬ 
tonic  structure  underneath  each  station  very  well,  mainly  due  to  the  broad  coverage  of 
azimuths  at  each  station.  Generally  speaking,  the  station  terms  are  positive  in  shield  regions 
such  as  Australia,  Canada,  India,  and  Scandinavia,  and  they  are  negative  in  the  east  Africa 
rift  valleys,  mid-ocean  ridges  {e.g. ,  Iceland  and  Azores  Islands),  island  arcs  {e.g. ,  Indonesia, 
Japan,  and  Taiwan),  and  Himalaya  Mountain  Ranges  (Chaman  Fault,  northern  India,  Nepal, 
and  Burman  Arc).  Solomon  and  Toksoz  (1970)  and  many  other  studies  {e.g.,  Evernden  and 
Clark,  1970;  Booth  etal.,  1974)  found  that  for  stations  in  U.S.,  the  attenuation  is  higher 
between  the  Rockies  and  Cascades,  and  In  the  northeastern  U.S.  This  pattern  is  also 
observable  in  Figure  1  (see  also  North,  1977).  As  North  (1977)  put  it,  it  is  gratifying  that  a 
simple  parameter  such  as  can  be  utilized  to  reveal  the  tectonics.  It  should  be  noted,  how¬ 
ever,  that  our  empirical  station  terms  also  include  the  effect  due  to  the  crustal  amplification  if 
such  local  site  effect  is  shared  by  all  ray  paths  from  different  test  sites  to  a  particular  station. 
This  could  be  the  reason  of  a  few  outliers  such  as  HNR  (Honiara,  Solomon  Islands),  PMG 
(Port  Moresby,  East  Papua  New  Guinea),  RAB  (Rabaul,  New  Britain),  and  BAG  (Baguio  City, 
Luzon,  Philippines)  which  do  not  show  negative  station  terms  as  would  be  expected  from  the 
strong  seismicity  in  that  region  {cf.  Figure  1).  Another  possible  reason  is  that  these  stations 
have  relatively  poorer  azimuthal  sampling  in  our  data  set,  and  hence  the  station  bias  at  these 
three  stations  is  not  well  constrained.  The  minor  discrepancy  between  the  deterministic 
corrections  by  Marshall  etal.  (1979)  and  our  empirical  corrections  could  be  due  to  the  same 
reason. 

Figure  2  shows  the  empirical  station  terms  determined  by  North  (1977)  and  Ringdal 
(1986)  with  ISC  recordings  of  world-wide  earthquakes  occurred  during  the  periods  of  1964- 
1973  and  1971-1980,  respectively.  After  the  mean  station  residual  of  0.104  m.u.  is  removed 
from  Ringdal's  station  terms,  the  spatial  pattern  of  his  corrections  appears  to  be  in  very  good 
agreement  with  that  of  North  (1977)  and  the  present  study. 

Figures  3  through  7  show  the  map  of  the  “pure  propagation  effect”  (top)  and  the  com¬ 
bined  station  amplification  (bottom)  defined  as  the  sum  of  the  receiver  term  and  the  path 
effect  for  explosions  detonated  in  five  source  regions  in  Eastern  Kazakhstan  which  include 
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Degelen  Mountain  [Deg]  and  Murzhik  [Mzk]  in  addition  to  the  three  subregions  defined  in 
Ringdal  and  Marshall  (1989);  southwestern  Balapan  [BSW],  northeastern  Balapan  [BNE],  and 
the  transition  zone  [BTZ]  between  BSW  and  BNE.  The  path  term  at  each  station  can  be 
regarded  as  the  azimuthal  variation  (towards  the  various  source  regions)  relative  to  the  aver¬ 
aged  station  amplification.  An  important  observation  is  that  all  these  five  test  sites  exhibit 
very  different  azimuthal  and  radial  amplitude  variations.  Degelen  and  Murzhik  events  are  sys¬ 
tematically  enhanced  in  the  western  U.S.  and  reduced  in  eastern  U.S.,  whereas  Balapan 
events  are  all  reduced  in  the  whole  U.S.  Murzhik  events  are  reduced  in  Scandinavia,  but 
Balapan  and  Degelen  events  get  enhanced  tfiere.  Such  highly  direction-dependent, 
distance-dependent,  and  site-dependent  patterns  of  the  amplitude  fluctuation  could  be  a 
diagnostic  for  the  path  effects  in  the  proximity  of  the  test  sites.  Back  projections  (e.g., 
Lynnes  and  Lay,  1990)  of  the  residuals  onto  the  upper  mantle  and  the  lower  crust  reveal 
that  similar  residuals  come  into  alignment  in  several  regions  partitioned  by  known  geologi¬ 
cal  features  (Jih  and  Wagner,  1991a).  Murzhik  events  recorded  in  the  western  U.S.  and  in 
northeast  Asia.  Degelen  events  in  the  western  U.S.,  and  SW  Balapan  events  at  western 
European  stations  must  pass  through  the  area  between  Chinrau  fault  and  Chingiz-Kalba 
shear  zone.  All  these  paths  show  positive  residuals.  The  area  north  of  Chinrau  fault 
might  have  some  complex  features  that  result  in  negative  mean  residuals.  Paths  from  NE 
Balapan  to  North  America  and  many  continental  European  stations  must  cross  this  area  or 
even  travel  along  the  Chinrau  fault  before  entering  the  deeper  mantle,  and  hence  the  com¬ 
plexity  in  the  waveforms  is  inevitable.  It  seems  that  the  mean  m^-Lg  separation  of  0.11-0.17 
m.u.  {e.g.,  Ringdal  and  Hokland,  1987;  Ringdal  and  Marshall,  1989;  Richards  etal.,  1990; 
Jih  and  Wagner,  1991a)  between  the  NE  and  SW  subregions  of  Balapan  could  be  due  in 
part  to  the  path  effects  —  in  addition  to  the  difference  of  source  medium  postulated  previously 
by  Marshall  et  al.  (1984).  A  detailed  discussion  on  the  seismic  variability  within  Balapan  test 
site  is  given  in  a  later  section  {cf.  pages  48-49).  Path  effects  can  also  explain  why  the  SW 
Balapan  waveforms  tend  to  be  more  complex  at  YKA  than  those  recorded  at  WRA,  EKA,  and 
GBA  arrays  (Jih  and  Wagner,  1991a). 

The  initial  P  waves  from  the  three  adjacent  test  sites  have  virtually  the  same  incident 
angle  at  each  teleseismic  station,  and  anything  in  common  across  all  events  (such  as  the  cru¬ 
stal  amplification  as  well  as  the  upper  mantle  attenuation  underneath  the  receiver)  would  have 
been  lumped  into  the  constant  station  term.  Thus  the  station  residuals  averaged  over  all 
events  from  the  same  test  site  would  correlate  very  little  with  the  receiver.  Instead,  they 
should  reveal  more  site-dependent  information  about  the  focusing/defocusing  pattern  under¬ 
neath  E.  Kazakhstan. 
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The  largest  and  most  prominent  fault  in  the  region  is  the  southeast-trending  Chingiz 
right-lateral  strike-slip  fault  that  passes  about  1 0  km  southwest  of  Degelen  Mountain  and  right 
across  the  Murzhik  test  area  (Rodean,  1979;  Bonham  eta!.,  1980;  Leith,  1987b).  Soviets 
reported  that  this  fault  has  a  very  steep  dip,  which  is  consistent  with  its  linear  expression  over 
large  distance  as  seen  on  Landsat  imagery  (Bonham  etai,  1980).  A  distinct  fault-line  scarp  is 
developed  along  much  of  the  oldest  metamorphic  rocks.  Chingiz  Fault  extends  for  a  total 
length  of  about  700  km.  Soviet  reports  postulate  that  this  fault  extends  down  to  the  boundary 
of  the  granite  layer  of  the  crust  and  possibly  into  the  upper  mantle.  For  Murzhik  explosions, 
the  propagation  of  P„  and  Lg  waves  could  be  affected  by  this  fault  significantly,  which  results 
in  a  radiation  pattern  such  as  we  observe  in  Figure  7.  More  specifically,  the  rays  towards  NW 
direction  could  be  reflected  or  diffracted  to  other  quadrants,  due  to  its  post-critical  incidence 
angles.  Such  relatively  distant  crustal  structure  should  have  little  impact  on  the  first  P  waves 
of  Balapan  explosions  at  teleseismic  distances,  however.  As  a  result,  amplitudes  of  Balapan 
events  recorded  at  Scandinavian  stations  are  still  largely  controlled  by  the  weak-attenuating 
shield  paths. 

Table  2  lists  the  WWSSN  station  corrections  and  path  corrections  for  explosions  in  the 
five  Central  Asian  nuclear  test  sites.  Note  that  the  station  terms  are  applicable  to  other  source 
regions  of  the  world  as  well. 
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EMPIRICAL  STATION  TERMS  VS.  DETERMINISTIC  STATION  TERMS 


Figure  1.  The  WWSSN  station  terms  (top)  inferred  from  a  GLM/MLE  joint  inversion  scheme  which 
simultaneously  inverts  for  the  seismic  source  sizes,  receiver  terms,  as  well  as  the  path  effects.  The 
inversion  of  2733  unknown  parameters  is  carried  out  with  13840  signals,  9080  noise  measurements,  and 
1609  clips  from  217  worldwide  explosions  recorded  at  118  selected  WWSSN  stations.  Only  paths  within 
20  and  95  degrees  are  used.  The  high  correlation  between  the  tectonic  type  and  the  station  terms  sug¬ 
gests  that  these  empirical  corrections  do  reflea  the  upper  mantle  conditions  underneath  the  receivers. 
Darkened  stars  represent  some  of  the  nuclear  test  sites  used  in  this  study.  The  match  between  the 
deterministic  station  corrections  derived  by  Marshall  et  al.  (1979)  (bottom)  and  our  empirical  corrections 
is  fairly  good.  14 
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EMPIRICAL  STATION  TERMS  INFERRED  FROM  ISC  EARTHQUAKE  RECORDINGS 


Figure  2.  The  empirical  station  terms  determined  by  North  (1977)  and  Ringdal  (1986)  with  world-wide 
earthquakes  recorded  at  ISC  during  the  periods  of  1964-1973  and  1971-1980,  respectively.  After  the 
mean  station  residual  of  0.104  m.u.  is  removed  from  Ringdal’s  station  terms,  the  spatial  pattern  of  his 
corrections  appears  to  be  in  very  good  agreement  with  that  of  North  (1977)  and  the  present  study. 
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Figure  3.  The  map  showing  the  “pure  propagation  effect”  (top)  and  the  combined  station  amplification 
(bottom)  defined  as  the  sum  of  the  receiver  term  (Figure  1)  and  the  path  effect  for  southwestern 
Balapan  explosions. 
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Figure  5.  Same  as  Figure  3  except  for  the  northeastern  region  of  Balapan  test  site. 
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Figure  6.  Same  as  Figure  3  except  for  Degelen  test  site. 


Table  2.  Receiver  and  Path  Effect  on  m/,  of  Semipalatinsk  Events 


Station 

Path  Terms  [F] 

Code 

Lon 

Lat 

Rev' 

BSW 

BTZ 

BNE 

Deg 

Mzk 

AAE 

38.766 

9.029 

-0.290 

-0.439 

-0.274 

-0.640 

-0.287 

-0.447 

AAM 

-83.656 

42.300 

0.210 

0.267 

0.140 

0.163 

-0.224 

-0.335 

AKU 

-18.107 

65.687 

-0.053 

0.169 

0.342 

0.155 

0.185 

0.043 

ANP 

121.517 

25.183 

-0.139 

-0.232 

-0.514 

0.123 

-0.301 

AQU 

13.403 

42.354 

-0.102 

-0.167 

-0.141 

-0.159 

-0.051 

-0.104 

ATU 

23.717 

37.972 

0.171 

0.050 

.\088 

0.040 

-0.312 

-0.176 

BAG 

120.580 

16.411 

0.030 

-0.189 

-0.148 

-0.064 

-0.213 

-0.273 

BEC 

-64.681 

32.379 

-0.120 

0.160 

0.216 

-0.072 

-0.087 

-0.289 

BKS 

-122.235 

37.877 

0.104 

-0.008 

-0.015 

-0.116 

-0.026 

0.112 

BLA 

-80.421 

37.21 1 

0.022 

-0.165 

-0.158 

-0.184 

-0.217 

-0.410 

B02 

-111.633 

45.600 

0.046 

-0.068 

-0.166 

0.055 

-0.074 

BUL 

28.613 

-20.143 

-0.004 

-0.112 

0.047 

-0.024 

-0.312 

-0.079 

CHG 

98.977 

18.790 

-0.234 

0.268 

0.623 

0.319 

0.177 

0.097 

CMC 

-115.083 

67.833 

-0.270 

0.134 

0.512 

0.522 

COL 

-147.793 

64.900 

-0.002 

0.041 

0.235 

0.220 

0.074 

-0.067 

COP 

12.433 

55.683 

0.174 

0.016 

0.032 

0.038 

0.047 

-0.499 

COR 

-123.303 

44.586 

0.161 

0.030 

0.127 

0.096 

0.095 

0.086 

CTA 

146.254 

-20.088 

0.130 

-0.096 

-0.026 

-0.074 

-0.081 

-0.130 

DAG 

-18.770 

76.770 

-0.076 

0.218 

0.100 

0.047 

0.169 

DAV 

125.575 

7.088 

-0.040 

-0.457 

-0.815 

-0.377 

-0.394 

DUG 

-112.813 

40.195 

0.074 

0.220 

0.056 

0.159 

0.390 

0.350 

EIL 

34.950 

29.550 

-0.067 

0.040 

-0.176 

-0.026 

ESK 

-3.205 

55.317 

0.084 

0.002 

-0.095 

-0.240 

0.085 

-0.522 

FLO 

-90.370 

38.802 

-0.100 

-0.110 

-0.005 

-0.518 

FVM 

-90.426 

37.984 

0.069 

-0.023 

-0.068 

-0.004 

-0.021 

GDH 

-53.533 

69.250 

-0.155 

0.180 

0.115 

0.103 

0.026 

0.235 

GEO 

-77.067 

38.900 

-0.006 

0.070 

0.024 

0.345 

-0.043 

-0.177 

1 )  mean  receiver  bias  that  needs  be  corrected  in  addition  to  the  path  correction. 
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Table  2.  Receiver  and  Path  Effect  on 


Station 


Code 

Lon 

Lat 

Rev 

GOL 

-105.371 

39.700 

-0.237 

GSC 

-116.805 

35.302 

0.022 

GUA 

144.912 

13.538 

-0.232 

HKC 

114.172 

22.304 

-0.087 

HLW 

31.342 

29.858 

-0.256 

HNR 

159.947 

-9.432 

0.220 

1ST 

28.996 

41.046 

0.184 

JER 

35.197 

31.772 

-0.014 

KBS 

11.924 

78.918 

-0.227 

KEV 

27.007 

69.755 

-0.108 

KOD 

77.467 

10.233 

0.177 

KON 

9.598 

59.649 

0.000 

KRK 

30.062 

69.724 

0.069 

KTG 

-21.983 

70.417 

-0.254 

LEM 

-6.833 

-0.417 

LON 

-121.810 

46.750 

-0.121 

LOR 

3.851 

47,267 

0.008 

MAL 

-4.41 1 

36.728 

-0.056 

MAN 

121.077 

14.662 

0.367 

MAT 

138.207 

36.542 

-0.256 

MDS 

-89.760 

43.372 

-0.091 

MSO 

-113.941 

46.829 

-0.091 

MUN 

116.208 

-31.978 

0.256 

NAI 

36.804 

-1 .274 

-0.110 

NDI 

77.217 

28.683 

0.100 

NHA 

109.212 

12.210 

-0.049 

NOR 

-16.683 

81.600 

-0.247 

BSW 


0.164 


0.058 


0.069 


0.015 


-0.077 


-0.197 


0.122 


-0.060 


-0.1 1 1 


0.320 


0.280 


0.352 


of  Semipalatinsk  Events 


Path  Terms  [F] 


BTZ  BNE  Deg 


0.187  0.144  0.181 


0.095  0.011  0.106 


-0.613  0.278 


-0.069  -0.017  -0.347 


0.006  0.004  -0.31 1 


0.111 


-0.059 


-0.341 


0.239 


0.119 


0.460 


-0.045 


0.099 


-0.361 


0.176 


0.213 


0.262 


0.132 


0.046 


0.002 


-0.244 


0.174 


-0.536 


-0.188 


0.002 


-0.019 


-0.184 


0.241 


-0.045 


0.032 


-0.138 


0.086 


0.098 


-0.528 


-0.188 


0.026 


0.137 


-0.097 


-0.188 


0.025 


-0.295 


-0.041 


0.175 


-0.279 


0.014 


0.013 


-0.004 


0.001 


-0.255 


-0.148 


0.058 


0.064 


0.045 


0.105 


0.030 


0.099 


-0.264 


0.143 


-0.122 


-0.020 


-0.220 


-0.166 


0.371 


-0.087 


-0.108 


-0.167 


0.008 


Mzk 


-0.050 


-0.104 


-0.179 


-0.235 


-0.015 


-0.457 


-0.064 


0.364 


-0.260 


0.044 


0.055 


0.067 


0.064 


-0.193 


-0.107 


-0.137 


-0.168 


0.134 


Table  2.  Receiver  and  Path  Effect  on  mb 


Station 


Code 


NUR 


OGD 


PDA 


PMG 


POO 


PRE 


PTO 


QUE 


RAB 


RCD 


SCP 


SDB 


SEO 


SHI 


SHK 


SHL 


SLR 


Lon 


24.651 


-74.596 


-25.663 


147.154 


73.850 


28.190 


-8.602 


66.950 


152.170 


-103.208 


-77.865 


13.572 


126.967 


52.520 


132.678 


91.883 


28.282 


Lat 


60.509 


41.088 


37.747 


-9.409 


18.533 


-25.753 


41.139 


30.188 


-4.191 


44.075 


40.795 


-14.926 


37.567 


29.638 


34.532 


25.567 


-25.735 


Rev 


0.089 


-0.164 


0.043 


0.177 


-0.004 


-0.074 


-0.193 


-0.484 


0.184 


0.334 


-0.060 


-0.049 


-0.125 


0.120 


-0.250 


-0.081 


-0.185 


BSW 


0.421 


-0.068 


-0.065 


-0.036 


0.126 


0.033 


-0.017 


0.090 


-0.258 


-0.016 


0.009 


-0.175 


-0.033 


-0.079 


0.169 


-0.268 


of  Semipalatinsk  Events 


Path  Terms  [F] 


BTZ  BNE  Deg  Mzk 


0.557  0.564  -0.036  0.006 


0.060  0.024  -0.258  -0.282 


0.005  -0.048  -0.316 


0.030  0.058  -0.016 


0.097  -0.016  0.131 


0.103  0.051  -0.232 


-0.081  -0.149  -0.121 


0.029  0.182  -0.095 


-0.160  -0.137  -0.591 


-0.259  -0.066  -0.232 


0.058  0.055  -0.131 


0.290  0.077  0.100 


-0.167  0.025  -0.452 


-0.038  0.083 


0.135  -0.377 


0.212  0.019  -0.056 


-0.343 


0.150 


0.266 


-0.100 


-0.084 


-0.226 


-0.341 


-0.264 


0.151 


-0.279 


0.002 


-0.514 


-0.117 


SNG 

100.620 

7.173 

0.005 

0.066 

0.109 

0.046 

-0.063 

0.008 

STU 

9.195 

48.772 

-0.001 

0.030 

0.049 

-0.024 

0.227 

0.153 

TAB 

46.327 

38.068 

0.290 

-0.074 

0.085 

0.170 

-0.003 

0.067 

TOL 

-4.049 

39.881 

0.120 

-0.072 

0.066 

-0.139 

-0.108 

-0.095 

TRI 

13.764 

45.709 

-0.193 

0.153 

0.134 

0.074 

0.271 

0.028 

UME 

20.237 

63.815 

0.070 

0.634 

0.302 

0.454 

-0.011 

-0.007 

VAL 

-10.244 

51.939 

-0.024 

0.050 

0.002 

0.105 

-0.105 

WES 

-71.322 

_ 

42.385 

-0.228 

-0.058 

0.013 

-0.352 

-0.357 
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Marshall  etai  (1992)  analyze  Degelen  and  Murzhik  events  recorded  at  4  U.K.-designed 
arrays,  and  they  find  that  EKA  and  GBA  have  distinguishable  path  effects  for  these  two  test 
sites.  Amplitudes  of  Murzhik  events  are  significantly  reduced  at  EKA,  whereas  those  of 
Degelen  events  are  magnified.  On  the  other  hand,  GBA  shows  a  strong  enhancement  for 
Muzhik  signals,  but  nearly  no  effect  on  Degelen  events.  At  YKA  or  WRA,  the  station/path 
effects  are  about  the  same  for  Degelen  and  Murzhik  explosions.  All  these  observations  (Fig¬ 
ure  6  of  Marshall  etai,  1992)  are  in  excellent  agreement  with  our  result  based  on  WWSSN 
recordings.  The  following  is  excerpted  from  Table  2,  which  illustrates  the  distinct  path  effects 
at  EKA  and  GBA  (see  also  Table  1  of  Jih  and  Wagner,  1991b): 


*)  A:  distance  from  EKA. 


Table  4.  Path  Terms  for  Stations  Close  to  GBA 

Station 

Test  Site 

« 

A 

Code 

Degelen 

Murzhik 

(km) 

KOD 

0.045 

0.364 

373 

POO 

0.131 

0.266 

666 

NDI 

0.008 

0.134 

1669 

*)  A;  distance  from  GBA. 


Note  that  the  consistent  trend  across  stations  of  wide  spatiai  spread  as  illustrated  in 
Tabies  3  and  4  suggests  that  these  path  effects  are  due  to  some  very  near-source 
focusing/defoucusing  feature. 


24 


To  further  illustrate  the  robustness  of  the  proposed  determination  scheme,  Figures  8 
through  15  show  the  raw  and  corrected  of  8  important  explosions  from  Eastern  Kaza¬ 
khstan.  Figure  8  shows  the  scatter  plot  of  3  different  types  of  station  m^s  for  Soviet  Joint 
Verification  Experiment  [JVE]  explosion  8809 14B.  The  25  good  recordings  and  1  clipped  sig¬ 
nal  are  shown  with  filled  circles  and  upward  arrow,  respectively.  The  raw  station  m^s  (top) 
have  a  standard  deviation  of  0.19  m.u.  Applying  the  “primary"  station  corrections  (/.e.,  the 
“Rev"  column  in  Table  2  or  the  “S"  term  in  (3))  and  the  “secondary"  corrections  (/.e.,  the 
“Path"  column  in  Table  2  or  “F"  term  in  [3])  reduces  the  scatter  down  to  0.11  m.u.  The 
dashed  lines  of  1  a  range  around  the  network-averaged  mt,  clearly  illustrate  the  remarkable 
reduction  of  fluctuation  across  the  recording  stations.  The  mean  event  itself  is  not 
significantly  changed,  however. 

Among  the  eight  Semipalatinsk  events  shown  here,  the  event  790804B  (Figure  9)  has 
the  smallest  scatter  in  the  resulting  m2.g  values.  The  dramatic  reduction  of  variation  from 
to  /772  9  shows  a  factor  of  nearly  3,  as  compared  to  the  worst  case  of  about  2.11  for  the  event 
811018B  (Figure  12).  Novaya  Zemlya  events  typically  exhibit  a  reduction  factor  around  2  (Jih 
and  Wagner,  1992).  Note  that  the  path  correction  proposed  in  this  study  not  only  reduces  the 

scatter  at  stations  that  reported  the  good  signals,  but  it  also  improves  the  data  con¬ 
sistency  of  the  censored  recordings,  as  indicated  by  the  shifting  of  the  clipping  recordings  In 
Figure  9. 
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VARIOUS  WWSSN  MAGNITUDES  OF  EVENT  88091 4B 


Azimuth  (degrees) 


Figure  8.  Scatter  plot  of  3  different  types  of  station  m^s  for  Soviet  JVE  {880914B).  The  25  good 
recordings  and  1  clip  are  shown  with  filled  circles  and  upward  arrow,  respectively.  The  raw  station  mj,s 
(top)  have  a  standard  deviation  of  0.19  m.u.  Applying  the  “primary"  station  corrections  and  the  “secon¬ 
dary"  path  corrections  reduces  the  scatter  down  to  0.11  m.u.  (txittom).  The  dashed  lines  around  the 
network-averaged  clearly  illustrate  the  remarkable  reduction  of  fluctuation  across  the  recording  sta¬ 
tions.  The  mean  event  rrit,  itself  is  not  significantly  changed,  however. 
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VARIOUS  WWSSN  MAGNITUDES  OF  EVENT  790804B 


0  SO  100  150  200  250  300  350 


0  50  100  150  200  250  300  350 


Azimuth  (degrees) 


Figure  9.  Scatter  plot  of  3  different  types  of  station  m^s  for  the  “historical"  Balapan  explosion  790804B. 
The  three  noise  measurements  are  shown  as  downward  arrows.  This  event  has  extremely  small  scatter 
in  the  resulting  m2g  values.  The  dramatic  reduction  of  variation  from  m,  to  m2g  shows  a  factor  of 
nearly  2.97. 


VARIOUS  WWSSN  MAGNITUDES  OF  EVENT  791028B 


Figure  10.  Scatter  plot  of  3  different  types  of  station  m/,s  for  the  “historical”  Balapan  explosion 
791028B. 
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VARIOUS  WWSSN  MAGNITUDES  OF  EVENT  791223B 


0  50  100  150  200  250  300  350 


0  50  100  150  200  250  300  350 


0  50  100  150  200  250  300  350 


A  zimuth  ( degrees) 


Figure  11.  Scatter  plot  of  3  different  types  of  station  robS  for  the  “historical"  Balapan  explosion 
7912238. 
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VARIOUS  WWSSN  MAGNITUDES  OF  EVENT  8110188 


Figure  12.  Scatter  plot  of  3  different  types  of  station  m^s  for  the  “historical"  Balapan  explosion 
811018B. 
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VARIOUS  WWSSN  MAGNITUDES  OF  EVENT  840526B 


Figure  13.  Scatter  plot  of  3  different  types  of  station  m^s  for  the  “historical"  Balapan  explosion 
840526B. 
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VARIOUS  WWSSN  MAGNITUDES  OF  EVENT  710425D 


Figure  14.  Scatter  plot  of  3  different  types  of  station  rObS  for  Degelen  explosion  710425D. 
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VARIOUS  WWSSN  MAGNITUDES  OF  EVENT  69122aM 


Figure  15.  Scatter  plot  of  3  different  types  of  station  m^s  for  Murzhik  explosions  691223M. 
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4.  YIELD  ESTIMATES  OF  SEMIPALATINSK  EXPLOSIONS 


It  is  fortuitous  to  have  the  source  information  released  by  Bocharov  etal.  (1989)  (and 
Vergino,  1 989)  to  calibrate  the  Semipalatinsk  test  site.  The  small  scatter  around  the  following 
calibration  curves  based  on  the  regression  of  our  path-corrected  on  the  published  yields 
illustrates  how  well  the  fit  can  be  at  the  Central  Asian  test  site. 

mbiPa)  =  0.802(±0.020)  log(W)  +  3.834(±0.032)  [6] 

mbiPb)  =  0.800(10.020)  log(W)  +  4.130(10.031 )  [7] 

mbiPmax)  =  0.768(10.019)  log(W)  +  4.399(10.030)  [8] 

Figures  16  through  18  show  the  regression  of  m2.g(Pa),  ^2.g{Pb)>  and  TT^2.Q(Pmax)  on  the 
the  Soviet  yields  published  by  Bocharov  efa/.  (1989),  which  correspond  to  Equations  [6] 
through  [8],  respectively.  The  uncertainties  in  the  m^s  and  the  yields  are  taken  into  account 
through  800  bootstrap  resamplings.  The  darkened  bundle  is  actually  the  collection  of  all  800 
regressions,  each  produced  by  a  possible  realization  of  11  perturbed  (m^,  yield)  pairs.  The 
95%  confidence  band  (shown  as  2  curves  around  the  darkened  bundle)  is  narrower  near  the 
centroid  and  wider  towards  both  ends,  as  expected.  The  individual  95%  confidence  intervals 
of  the  two  inferred  parameters  (/.e. ,  the  slope  and  the  intercept  of  the  calibration  curve)  are 
shown  with  the  dashed  line  in  the  scatter  plot  (bottom).  Note  that  the  dashed  rectangle  is  not 
the  joint  90%  confidence  interval,  however,  due  to  the  highly  correlated  nature  of  the  two 
parameters.  Degelen  event  660507D  is  not  included  in  these  regressions,  as  suggested  by 
Jih  and  Wagner  (1991b). 

We  have  utilized  these  calibration  curves  to  estimate  the  yield  of  all  92  Semipalatinsk 
explosions  in  our  data  set,  and  the  result  is  summarized  in  Table  5.  For  cratering  events 
(such  as  6501 15B)  the  yield  estimate  based  on  the  first  motion  (/.e. ,  Pg)  should  be  used, 
since  no  depth  correction  (e.g. ,  Marshall  efa/.,  1979)  has  been  applied  to  mb{Pb)  or 
iTibiPmax)  '0  Table  5.  For  this  particular  event,  Myasnikov  etal.  (1970)  gave  a  "scaled 
apparent  radius"  and  scaled  depth  of  of  51  and  50  m/KT°^^,  respectively.  Combining  this 
information  with  the  crater  radius  and  the  emplacement  depth  released  at  the  IAEA  sympo¬ 
sium,  Ringdal  and  Marshall  (1989)  inferred  the  yield  of  this  explosion  as  111  KT,  which  is 
identical  to  our  estimate  based  on  Pg  (Table  5).  This  example  illustrates  that  Pg  from  hard- 
rock  test  sites  in  stable  region  could  be  a  very  favorable  phase  for  the  source  size  determina¬ 
tion. 
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Much  of  the  source  information  about  the  Soviet  JVE  explosion  (8809148,  Figure  8)  has 
not  been  released.  The  "New  York  Times"  (Gordan,  1988)  states  that  the  American  and 
Soviet  on-site  measurements  are  said  to  give  yields  of  115  KT  and  122  KT,  respectively.  If 
we  substitute  the  Tn2Q{Pa)  of  JVE  into  Equation  [6],  the  mean  yield  estimate  would  be  1 13  KT 
with  a  95%-uncertainty  factor  of  2.2  {cf.  Figure  16).  If  the  averaged  on-site  measurement  of 
119  KT  was  accurate,  then  our  central-value  estimate  of  113  KT  based  on  P-wave  observa¬ 
tions  alone  is  off  by  5%  only. 

The  Soviet  JVE  has  a  RMS  Lg  measured  at  NORSAR  of  5.969±0.01  (Ringdal  and 
Marshall,  1989).  Along  with  the  other  8  RMS  Lg  values  in  Ringdal  (1990),  DWLSQ-derived 
calibration  curve  would  give  an  estimate  of  1 03  KT  for  JVE  with  a  95%-factor  of  1 .8.,  assum¬ 
ing  a  10%  standard  error  in  tfie  yields  (Figure  14  of  Jih  and  Wagner,  1991b).  This  yield  esti¬ 
mate  is  nearly  Identical  to  the  one  obtained  with  the  in-country  regional  network  recordings. 
The  RMS  Lg  of  5.980±0.02  (Israelson,  1991)  corresponds  to  100  KT  and  a  95%-factor  of  1.9 
based  on  either  DWLSQ  of  Jih  et  al.  (1991)  or  Ericsson’s  algorithm.  In  addition  to  the  yield 
estimates  cited  above,  Sykes  and  Ekstrom  (1989)  gave  an  estimate  of  113  KT  based  on  the 
arithmetic  average  of  m^,  and  Mg.  Priestley  et  aL  (1990)  analyzed  the  Lg  amplitudes  at  4 
seismographs  near  the  Semipalatinsk  test  range;  KSU  (Karasu),  KKL  (Karkaralinsk),  BAY 
(Bayanual),  and  TLG  (Talgar),  and  they  obtained  a  m^iLg)  of  5.968±0.02.  Murphy  et  al. 
(1991)  gave  a  network-averaged  /i?/,  of  6.012  with  a  standard  deviation  of  0.190  across  the 
network.  They  also  derived  a  RMS  Lg  of  5.969  using  8  stations  in  U.S.S.R.,  Norway,  and 
Manchuria.  It  is  worth  noting  that  all  these  seismic  magnitudes  give  very  consistent  yield  esti¬ 
mate  in  the  range  100-150  KT,  as  specified  in  the  bilateral  agreement  signed  by  U.S.  and 
Soviet  governments  before  JVE  (Richards,  1990;  Stump,  1991). 

There  are  15  events  in  common  in  Israelson’s  (1991)  RMS  Lg  and  our  7^2,9  data  sets  for 
which  the  Soviet-published  yields  are  available,  which  yield  a  very  weak  correlation  between 
the  RMS  Lg  and  nia  g  residuals  (relative  to  the  expected  magnitude  at  the  associated  yield 
value)  (Figure  19),  and  hence  the  combination  of  these  two  methods  for  better  yield  estimate 
is  justifiable. 
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mb:Yield  Relationship  at  KTS 


OWLS  (uncertain  X  &  Y):  8=0^(0.020),  1=3.83(0.032),  18.  data  used, 

95%  error  In  mb  at  1,1 0,50,1 00,1 50KT:  0.24, 0.13, 0.11, 0.13, 0.16, 

95%  factor  In  yield  at  1,1 0,50,1 00,1 50ia:  380, 2.06, 1.92, 2.17, 2.44 

OWLS  (precise  X  assumed):  8=083(0.043),  1=380(0.067) 

Standard  LS:  8=0.81(0.041),  1=3.83(0.064) 

10%  S.E.  In  yields  assumed 

Scatter  Plot  of  Inferred  Parameters 


95%  confidence  Interval  of  slope:  0.80W-0.043 
95%  confidence  Interval  of  Intercept;  3.83+/-0.068 
[97.5%  quantile  of  t(l  6.  D.o.F.),  2.120,  used] 


Figure  16.  Regressing  the  on  19  Soviet-published  yields.  The  yields  are  assumed  to  be  subject 

to  10%  standard  errors.  The  uncertainties  in  the  m^s  and  the  yields  are  taken  into  account  through  800 
Monte-Carlo  resamplings.  The  darkened  bundle  is  actually  the  collection  of  all  800  regressions,  each 
associated  with  a  possible  realization  of  19  perturbed  (mj,,  yield)  pairs.  The  95%  confidence  band 
(shown  as  2  curves  around  the  darkened  bundle)  is  most  narrow  near  the  centroid  and  wider  towards 
both  ends,  as  expected.  The  individual  95%  confidence  intervals  of  the  two  inferred  parameters  (i.e. ,  the 
slope  and  the  intercept  of  the  calibration  curve)  are  shown  with  the  dashed  line  in  the  scatter  plot  (bot¬ 
tom).  Note  that  the  dashed  rectangle  is  not  the  joint  90%  confidence  interval,  however,  due  to  the 
highly  correlated  nature  of  the  two  parameters. 


800-bootstrap  Intercepts  nib(WWSSN,MLE8,Pb) 


nnb:Yield  Relationship  at  KTS 


OWLS  (uncertain  X  &  Y):  8=0.80(0.020),  1=4.13(0.031),  19.  data  used, 

95%  error  In  irib  at  1,1 0,50,1 00,1 SOKT:  0.22, 0.12, 0.10, 0.12, 0.13, 

95%  factor  In  yield  at  1 ,1 0,50,1 00,1 SOKT:  3.45, 2.01 , 1 .81 , 2.02, 2.1 7 

OWLS  (precise  X  assumed):  8=0.82(0.035),  1=4.10(0.055) 

Standard  LS:  8=0.80(0.032),  1=4.12(0.052) 

10%  S.E.  In  yields  assumed 


Scatter  Plot  of  Inferred  Parameters 
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95%  confidence  Interval  of  slope:  0.804^/-0.041 
95%  confidence  Interval  of  intercept:  4.13+/-O.065 
(97.5%  quantile  of  t(1 7.  D.o.F.),  2.110,  used) 


Figure  17.  Same  as  Figure  16  except  for  mtJiPb) . 


800-bootstrap  Intercepts  mb(WWSSN,MLE8,Pmax) 


mb:Yield  Relationship  at  KTS 


OWLS  (uncertain  X  &  Y):  3=0.77(0.019),  1=4.40(0.030),  19.  data  used, 

95%  error  In  mb  at  1,1 0,50,1 00,1 50KT:  OJM,  0.12, 0.10, 0.12, 0.13, 

95%  factor  In  yield  at  1,1 0,50,1 00,1 50KT:  3.39, 2.01, 1.81, 2.02, 2.17 

OWLS  (precise  X  assumed):  3=0.78(0.032),  1=4.38(0.051) 

Standard  LS:  3=0.77(0.030),  1=4.39(0.048) 

10%  S.E.  In  yields  assumed 

Scatter  Plot  of  Inferred  Parameters 


95%  confidence  Interval  of  slope:  0.77+/-0.040 
95%  confidence  Interval  of  Intercept:  4.40+/-0.063 
(97.5%  quantile  of  t(l  7.  D.o.F.),  2.110,  used] 


Figure  18.  Same  as  Figure  16  except  for 
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Table  5.  Yield  Estimates  of  Semipalatinsk  Explosions 


Event 


Epicenter 


Yield  Estimate 


6501 15B 

BTZ 

79.009 

49.935 

111 

94 

81 

6511 21 D 

Deg 

78.064 

49.819 

25 

24 

24 

66021 3D 

Deg 

78.121 

49.809 

223 

197 

192 

660320D 

Deg 

78.024 

49.762 

94 

91 

94 

660507D 

Deg 

78.105 

49.743 

2 

1 

1 

661019D 

Deg 

78.021 

49.747 

46 

41 

36 

661218M 

Mzk 

77.747 

49.925 

88 

75 

78 

670226D 

Deg 

78.082 

49.746 

100 

89 

94 

67091 6M 

Mzk 

77.728 

49.937 

11 

10 

10 

670922M 

Mzk 

77.691 

49.960 

7 

8 

9 

671122M 

Mzk 

77.687 

49.942 

1 

1 

68061 9B 

BNE 

78.986 

49.980 

11 

12 

13 

680929D 

Deg 

78.122 

49.812 

54 

53 

51 

690531 M 

Mzk 

77.694 

49.950 

6 

9 

9 

690723D 

Deg 

78.130 

49.816 

12 

13 

13 

69091  ID 

Deg 

77.997 

49.776 

2 

2 

3 

691130B 

BTZ 

78.956 

49.924 

80 

101 

94 

691228M 

Mzk 

77.714 

49.937 

61 

60 

58 

700721 M 

Mzk 

77.673 

49.952 

12 

13 

14 

701104M 

Mzk 

77.762 

49.989 

24 

18 

17 

710322D 

Deg 

78.109 

49.798 

40 

40 

35 

710425D 

Deg 

78.034 

49.769 

99 

91 

88 

710606M 

Mzk 

77.660 

49.975 

20 

23 

22 

Yield 


Announced 


100-150 


29 


125 


100 


20-150 


20-150 


20-150 


<20 


10 


<20 


<20 


60 


<20 


16 


<20 


125 


46 


<20 


<20 


20-150 


90 


16 


Table  5.  Yield  Estimates  of  Semipaiatinsk  Explosions  (continued) 


Event 


Epicenter 


Yield  Estimate 


Date 

Site 

Lon 

Lat 

710619M 

Mzk 

77.641 

49.969 

710630B 

BTZ 

78.980 

49.946 

711009M 

Mzk 

77.641 

49.978 

711 021 M 

Mzk 

77.597 

49.974 

711230D 

Deg 

78.037 

49.760 

72021 OB 

BNE 

78.878 

50.024 

720328D 

Deg 

78.076 

49.733 

72081 6D 

Deg 

78.059 

49.765 

720826M 

Mzk 

77.717 

49.982 

720902M 

Mzk 

77.641 

49.959 

721102B 

BSW 

78.817 

49.927 

721210D 

Deg 

78.058 

49.819 

721210B 

BNE 

78.996 

50.027 

730723B 

BTZ 

78.850 

49.980 

731214B 

BNE 

79.010 

50.040 

750427B 

BNE 

78.980 

49.990 

760704B 

BTZ 

78.950 

49.910 

761123B 

BNE 

79.000 

49.990 

761207B 

BSW 

78.900 

49.880 

770329D 

Deg 

78.140 

49.790 

770730D 

Deg 

78.160 

49.770 

780326D 

Deg 

78.070 

49.730 

780422D 

Deg 

78.170 

49.720 

Yield 


Announced 


<20 


<20 


12 


23 


20-150 


16 


6 


8 


<20 


2 


165 


20-150 


140 
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Table  5. 


Event 


Yield  Estimates  of  Semipalatinsk  Explosions  (continued) 


Epicenter 


Yield  Estimate 


Date 

Site 

Lon 

Lat 

780611 B 

BSW 

78.838 

49.879 

780705B 

BSW 

78.871 

49.887 

780728D 

Deg 

78.140 

49.756 

780829B 

BNE 

78.990 

50.000 

78091 5B 

BTZ 

78.940 

49.910 

781104B 

BNE 

78.943 

50.034 

781129B 

BSW 

78.760 

49.950 

790623B 

BTZ 

78.910 

49.910 

790707B 

BNE 

79.060 

50.050 

790804B 

BSW 

78.904 

49.894 

79081 8B 

BNE 

79.010 

49.970 

791028B 

BNE 

78.997 

49.973 

791202B 

BSW 

78.840 

49.890 

791223B 

BSW 

78.755 

49.916 

800522D 

Deg 

78.082 

49.784 

800629B 

BSW 

78.815 

49.939 

80091 4B 

BTZ 

78.880 

49.970 

801012B 

BNE 

79.080 

49.950 

801214B 

BTZ 

79.000 

49.930 

801227B 

BNE 

79.040 

50.040 

810422B 

BSW 

78.900 

49.900 

810913B 

BTZ 

78.980 

49.890 

811018B 

BSW 

78.859 

49.923 

Yield 


Announced 
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Table  5.  Yield  Estimates  of  Semipalatinsk  Explosions  (continued) 


Yield  Estimate 


1  Event 

Epicenter 

Date 

Site 

Lon 

Lat 

811129B 

— 

BSW 

78.860 

49.887 

811227B 

BSW 

78.870 

49.900 

820425B 

BTZ 

78.930 

49.880 

820704B 

BTZ 

78.850 

49.990 

820831 B 

BSW 

78.761 

49.924 

821205B 

BSW 

78.840 

49.910 

821226B 

BNE 

78.988 

50.071 

83061 2B 

BTZ 

78.980 

49.910 

831006B 

BSW 

78.840 

49.930 

831026B 

BTZ 

78.910 

49.920 

831120B 

BNE 

78.999 

50.047 

840425B 

BTZ 

78.940 

49.950 

840526B 

BNE 

79.006 

49.969 

8407 14B 

BTZ 

78.960 

49.890 

841027B 

BSW 

78.830 

49.950 

841202B 

BNE 

79.070 

49.990 

841216B 

BTZ 

78.860 

49.960 

841228B 

BSW 

78.750 

49.860 

85021 OB 

BSW 

78.781 

49.888 

85061 5B 

BSW 

78.880 

49.890 

850630B 

BSW 

78.658 

49.848 

870620B 

BSW 

78.740 

49.927 

88091 4B 

BSW 

78.808 

49.833 

Yield 


Announced 
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5.  MISCELLANEOUS  COMPARATIVE  STUDIES  WITH 


In  the  initial  attempt  of  Jih  and  Wagner  (1991b)  to  correct  for  the  path  effects  of  Semi- 
palatinsk  explosions  with  "all  three"  subregions  of  Balapan  test  site  treated  as  geologically 
and  geophysically  homogeneous,  they  reported  that  3  out  of  82  Semipalatinsk  explosions 
appear  anomalous  in  that  tti^  q  do  not  show  significant  reduction  in  the  fluctuational  variation 
across  the  recording  network.  The  Soviet  JVE  explosion  was  one  of  the  three  anomalous 
events.  It  turns  out  that  separating  Balapan  test  site  into  three  regions  as  suggested  in  Ring- 
dal  and  Marshall  (1989)  (and  Marshall  etai,  1984;  Ringdal  and  Fyen,  1988)  improves  the 
performance  significantly,  and  none  of  the  three  events  would  seem  anomalous.  Such  parti¬ 
tioning  of  a  test  site  according  to  the  geophysical  characteristics  within  a  test  site  can  provide 
very  accurate  yield  estimate,  as  evidenced  by  our  yield  estimation  exercise  using  Pa  and  Pt, 
{cf.  pages  34-35). 

Tables  6  and  7  compare  P^ax  and  relative  to  Pg  at  several  Central  Asian  nuclear  test 
sites.  Note  that  there  appears  to  be  a  bias  of  0.11  m.u.  in  mt,(P^ax) -mb(Pa)  between  Eastern 
Kazakh  and  Novaya  Zemlya.  This  bias  could  be  largely  due  to  the  difference  in  pP  interfer¬ 
ence  at  these  two  test  sites,  however  (Jih  and  Wagner,  1992). 


Table  6.  m2.9(ML8),  and  mt,(Pt,)  vs.  mt,{Pg) 

Site^ 

my^Pb)-mb{Pa) 

^b(^max)  ^biPa) 

#2 

BSW 

0.288±0.013 

0.52310.014 

13 

BTZ 

0.289±0.033 

0.51610.059 

4 

BNE 

0.236±0.020 

0.44510.013 

5 

Deg 

0.317±0.027 

0.55710.037 

6 

KTS 

0.285±0.010 

0.51510.013 

28 

NNZ 

0.224±0.012 

0.40210.014 

17 

PRC 

— 

0.099±0.051 

— 

0.29410.073 

9 

1)  BSW  =  SW  subsite,  Balapan;  BNE  =  NE  subsite,  Balapan;  B7Z  =  transition  zone,  Balapan;  Deg  =  Degelen 
Mountain;  Mzk  »  Murzhik;  KTS  -  all  5  subsites  in  Eastern  Kazakh  combined;  NNZ  =  northern  subsite,  Novaya 
Zemlya;  PRC  =  Lop  Nor,  Tarim  Basin. 

2)  only  events  after  04/01/76  are  used. 
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Table  7.  TG(ML8),  m^{P^)  and  m<,(Pfc)  vs.  m,,(P^ 

Site 

mb{Pb)-mt{Pa) 

^b{Pmax)  ~^b{Pa) 

BSW 

0.290±0.012 

0.52410.013 

14 

BTZ 

0.290±0.021 

0.503K3.035 

8 

BNE 

0.268±0.018 

0.49310.025 

9 

Deg 

0.28910.012 

0.51510.014 

21 

Mzk 

0.30010.017 

0.53110.020 

13 

KTS 

0.28910.006 

0.51610.008 

65 

NNZ 

0.22610.008 

0.39610.009 

28 

PRC 

0.13310.042 

0.32810.063 

12 

1 )  all  events  used. 


Before  Bocharov  etal.  (1989)  published  the  yields  and  other  source  information  of  his¬ 
torical  Soviet  events,  several  attempts  had  been  made  to  investigate  the  characteristics  of 
cratering  explosions  in  that  region.  For  instance,  McLaughlin  etal.  (1985)  studied  the  ratio  of 
the  Pa  phase  and  P^ax  phase  of  presumed  Balapan  contained  and  cratering  explosions  by 
comparing  the  WWSSN  station  m^’s.  The  motivation  was  that  the  logarithm  of  amplitude  ratio 
of  Pmax^Pa  0^  event  6501 15B  was  significantly  smaller  than  other  presumed  contained  explo¬ 
sions  in  the  vicinity.  Assuming  the  phase  Pg  is  unaffected  by  the  influence  of  the  non-linear 
free-surface  interference,  then  an  adjustment  to  the  mbiPmwt)  should  be  able  to  convert  that 
to  a  contained  explosion  of  the  same  yield.  McLaughlin  et  a!.  (1985)  concluded  that  a  correc¬ 
tion  between  0.17  and  0.27  is  needed  for  this  conversion,  assuming  a  yield  of  125KT. 

Der  etal.  (1985)  deconvolved  four  contained  and  the  cratering  Balapan  events  6501 15B 
recorded  at  EKA,  and  then  they  convolved  the  Green’s  functions  with  an  appropriate  attenua¬ 
tion  operator  as  well  as  the  source-time  function  of  various  yields  of  interest.  By  comparing 
the  phases  Pg  and  P^ax  of  the  synthetics,  they  obtained  a  cratering-to-contained  correction  of 
0.15,  0.15,  and  0.18  at  60,  125,  and  300KT,  respectively. 

Day  et  al.  (1986)  did  a  theoretical  study  with  nonlinear  source  calculations  to  account  for 
coupling  variations  with  depth.  Their  results  are  summarized  as  follows; 
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n7b  (contained)  -  mt,  (cratering)  =  -  (0.1  to  0  15)  direct  coupling 

+  (0.1  to  0.25)  surface  interaction  effects 
=  (0  to  0.15)  total  bias. 

Based  on  46  Balapan  explosions  recorded  at  EKA,  Ringdal  and  Marshall  (1989)  derived 
a  value  of  0.75  as  their  mean  \og(P^ax/Pg)  across  the  EKA  array  using  the  same  techniques 
as  used  in  McLaughlin  et  al.  (1985).  The  cratering  event  6501 15B  had  nibiPmax)  -  ^b(Pa)  = 
0.62  at  EKA,  and  hence  they  apply  a  correction  of  5.87  +  (0.75  -  0.62)  =  6.00  for  a  hypotheti¬ 
cal  contained  explosion  with  equivalent  yield. 

We  utilize  the  statistics  in  Tables  6  through  9  to  illustrate  that  the  correction  by  Ringdal 
and  Marshall  (1989)  might  be  slightly  more  accurate  than  that  in  the  other  studies  cited 
above.  For  event  6501 15B,  our  mi,(P^ax)  *  ^b(^a)  =  0.392  (Table  1),  which  is  0.123  m.u. 
lower  than  the  average  mt,{Pmay)  -  mij{Pg)  =  0.515  (KTS)  or  0.516  (BTZ,  where  this  event  is 
located)  shown  in  Table  6.  This  implies  a  corrected  mgg  for  contained  111  KT  as 
5.988±0.013.  Alternatively,  using  Equations  [6]  and  [8],  we  expect  the  mtjl^P^gx)  -  rribiP^  to 
be  0.497  and  0.490  m.u.  at  100  and  150  KT,  respectively.  This  would  imply  a  correction  of 
about  0.1  m.u.  associated  with  the  range  of  100-150  KT. 

The  cratering-to-contained  conversions  cited  above  typically  require  extra  information 
about  the  general  behavior  of  contained  explosions  in  the  same  source  region.  For  the  pur¬ 
pose  of  estimating  the  yield  of  a  cratering  shot  in  an  isolated  region,  using  Pg  could  be  a 
much  easier  approach. 

Table  8  lists  the  expected  mi,  values  for  each  of  Pg,  Pi,,  and  P^ax  phases  from  NTS, 
NNZ,  and  KTS  explosions  based  on  Equations  [6]  and  [8].  The  estimated  "mean"  bias 
can  then  be  computed  in  a  straightforward  manner  (Table  9).  The  bias  estimates  based  on 
7772,2  are  included  for  comparison. 
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Table  8.  Expected 


7712,2 


m.  at  Various  Sites 

b 


Phase/Site 

1KT 

10KT 

50KT 

100KT 

1KT 

10KT 

50KT 

100KT 

4.270 

5.080 

5.646 

5.890 

4.399 

5.167 

5.704 

5.935 

TTiiJiPb)  (KTS) 

4.018 

4.848 

5.428 

5.678 

4.130 

4.930 

5.489 

5.730 

mb{Pa){^S) 

3.735 

4.560 

5.137 

5.385 

3.834 

4.636 

5.197 

5.438 

mb(P,„ax){NNZ)' 

4.254 

5.019 

5.554 

5.784 

4.245 

5.040 

5.596 

5.835 

(NNZ)' 

3.948 

4.785 

5.371 

5.623 

3.989 

4.835 

5.426 

5.681 

mb(Pa)(NNZ)’ 

3.735 

4.560 

5.137 

5.385 

3.834 

4.636 

5.197 

5.438 

^b(Pmax)(NTS) 

3.749 

4.608 

5.208 

5.467 

3.954 

4.771 

5.342 

5.588 

mb{Pb)mS) 

3.479 

4.349 

4.957 

5.219 

3.674 

4.505 

5.086 

5.336 

mt,(Pa)(NTS) 

3.368 

4.186 

4.758 

5.004 

3.607 

4.372 

4.907 

5.137 

1 )  from  Jih  and  Wagner  (1 992). 


Table  9.  Expected  mj,  Bias  Relative  to  NTS 

7712,2 

7712,9 

Phase/Site 

10KT 

50KT 

100KT 

150KT 

1CKT 

50KT 

100KT 

150KT 

max)  (KTS) 

0.47 

0.44 

0.42 

0.41 

0.40 

0.36 

0.35 

0.34 

mfc{Pi,)(KTS) 

0.50 

0.47 

0.46 

0.45 

0.42 

0.40 

0.39 

0.39 

mb{P^{KTS) 

0.37 

0.38 

0.38 

0.38 

0.26 

0.29 

0.30 

0.31 

m6(P„ax)(NNZ)’ 

0.41 

0.35 

0.32 

0.30 

0.27 

0.25 

0.25 

0.24 

mb(Pb)(NNZ)^ 

0.44 

0.41 

0.40 

0.40 

0.33 

0.34 

0.35 

0.35 

m,{P,)  {mz)' 

0.37 

0.38 

0.38 

0.38 

0.26 

0.29 

0.30 

0.31 

1)  from  Jih  and  Wagner  (1992). 


Marshall  et  al.  (1984)  found  that  explosions  in  the  northeast  and  southwest  portions  of 
Balapan  test  site  produce  distinctly  different  waveforms  when  recorded  at  the  UK  seismologi- 
cal  array  stations,  suggesting  that  Balapan  test  site  can  be  subdivided  into  two  areas 
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characterized  by  different  geophysical  properties.  Ringdal  and  HokJand  (1987)  find  that  this 
pattern  is  persistently  present  whether  based  on  worldwide  network  or  ^^(Pcoda)  of 
NORSAR  is  used.  They  Inferred  the  average  mt,-Lg  between  SW  and  NE  subregions  as  0.17 
m.u.  In  a  follow-up  study,  Ringdal  and  Fyen  (1988)  suggest  that  there  appear  to  be  a  transi¬ 
tion  zone  between  the  NE  and  SW  subregions.  Ringdal  and  Marshall  (1989)  recomputed  the 
SW-NE  bias  as  0.15  m.u.  with  96  Balapan  events  recorded  at  ISC  stations  and  NORSAR. 
Although  Ringdal  and  Marshall  (1989)  agree  that  the  possibility  of  a  m^(L^)  bias  contributing 
to  this  difference  between  SW  and  NE  cannot  be  entirely  ruled  out,  they  tend  to  believe  that 
this  bias  Is  due  to  a  relative  /n^,  bias  between  these  two  areas. 

We  followed  the  zoning  in  Ringdal  and  Marshall  (1989)  in  partitioning  Balapan  test  site 
into  three  regions:  southwest  (SW),  transition  zone  (TZ),  and  northeast  (NE).  The  mt,-Lg 
value  of  0.11  m.u.  shown  in  Table  10  is  slightly  smaller  than  that  of  previous  studies,  suggest¬ 
ing  that  our  ^2.9  (^max)  niay  be  somewhat  better  than  other  mj,.  Regressions  with  yields  pub¬ 
lished  by  Bocharov  et  al.  (1989)  show  that  NE  explosions  have  positive  Lg  residuals  and 
negative  rrit,  residuals,  whereas  SW  explosions  show  the  opposite  trend  (Figure  19).  Thus  it 
would  seem  plausible  that  the  apparent  mt^-Lg  bias  could  have  been  “enhanced”  by  the 
negative  correlation  between  mt,  and  Lg  residuals.  It  is  interesting  to  note  the  much  smaller 
mt,-Lg  bias  when  Pg  is  used  (Table  10).  A  three-dimensional  geological  model  of  the 
Balapan  test  site  by  Leith  and  Unger  (1989)  shows  a  distinct  difference  between  the  NE  and 
SW  portions  of  the  test  site,  with  the  granites  closer  to  the  surface  and  the  alluvium  thinner  in 
the  southwest.  The  thicker  alluvium  layer  in  NE  region  could  increase  the  waveform  complex¬ 
ity  and  reduce  the  magnitudes  measured  with  P^ax-  Nevertheless,  the  first  motion  (/.e. ,  Pg) 
should  be  least  affected  by  this  factor,  and  therefore  a  favorable  source  measure  —  so  long 
as  it  is  not  contaminated  by  the  microseismic  noise  at  the  receiver  site. 

Nuttli  (1987,  1988)  suggests  that  there  is  a  bias  of  about  0.2  m.u.  between  Degelen 
and  Balapan,  with  Degelen  explosions  having  even  larger  rrit,  excitation  (relative  to  Lg).  We 
do  not  see  such  Degelen-Balapan  bias  with  RMS  Lg  measured  at  NORSAR  (Table  10).  The 
Degelen  data  set  alone  is  too  small  for  decisive  conclusion.  However,  if  we  treat  Murzhik  as 
part  of  Degelen,  as  did  Nuttli  (1987),  the  average  mt,(P^gx)-RMS  L^(NORSAR)  bias  between 
Degelen  and  Balapan  is  only  0.02  m.u.,  which  is  insignificant. 
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Table  10.  m2g-RMS  Lp(NORSAR)  at  Various  Sites 


Site 

mb{Pa)-mb{l-g),  # 

mb{Pb)-mb(Lg),  # 

^fc(^max)  1  ^ 

BSW 

-0.504±0.011  11 

-0.228±0.011  11 

-i-0.023l0.01 5  20 

BTZ 

-0.523±0.045  6 

-0.243i<).020  6 

-0.04110.015  14 

BNE 

-0.565±0.023  8 

-0.304±0.014  9 

-0.09210.012  14 

Deg 

-0.484±0.046  5 

-0.207±0.042  5 

-i-0.012l0.034  5 

Mzk 

-0.562±0.073  3 

-0.25910.045  3 

-0.04610.032  3 

KTS 

-0.524±0.013  33 

-0.25010.010  34 

-0.02610.010  56 

NNZ 

-0.51 9±0 .020  14 

-0.29610.023  14 

-0.12110.024  14 

Table  1 1 .  ^2.9  * 

miJi^Lg)  (Nuttli)  at  Various  Sites 

Site 

mb{Pa)-mb(Lg).  # 

mdPb)-nibiLg).  # 

^biPfuax)  ~^l^^g)  •  ^ 

BSW 

-0.54410.050  8 

-0.244l0.040  8 

-0.03410.028  16 

BTZ 

-0.44410.078  5 

-0.18410.057  5 

-0.03310.024  13 

BNE 

-0.58410.043  6 

-0.30310.039  7 

-0.07510.031  14 

Deg 

-0.52410.124  5 

-0.19610.1 12  5 

-h0.05010.100  5 

KTS 

-0.52910.033  24 

-0.23910.028  25 

-0.03710.017  48 

NNZ 

-0.54910.033  24 

-0.32410.035  24 

-0.1 51 10.032  24 

We  suggest  that  the  mt,-Lg  bias  between  SW  and  NE  Balapan  can  be  tentatively 
decomposed  into  several  parts; 

[I]  Difference  in  pP  between  SW  and  NE, 

[II]  Difference  in  mi,  coupling,  i.e.,  mt,{S\N)  >  r77(,(NE), 

[III]  Difference  in  Lg  coupling,  i.e.,  Lg(HE)  >  Lg(SW), 

[IV]  Effects  due  to  the  station -station  correlation  structure, 

[V]  Effects  due  to  the  uneven  geographical  clustering  of  stations,  as  well  as  any  path  effect 
which  is  not  fully  accounted  for  through  the  network  averaging. 
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Based  on  our  [I]  is  about  0.05  m.u.,  whereas  [II]  and  [III]  are  about  0.02-0.03  m.u. 
each.  The  bias  of  0.11  m.u.  in  Table  10  using  mt,{Pmax)  'S  the  sum  of  [I]  through  [III].  It 
reduces  to  0.06  if  based  on  the  first  motion  is  used,  which  reflects  the  sum  of  [II]  and  [III]. 
For  ISC  data,  we  estimate  that  [V]  is  about  0.02  m.u.  if  m22  derived  by  the  conventional 
LSMF  are  used.  When  g  is  used,  this  term  is  eliminated,  and  hence  a  smaller  mi,-Lg  bias 
is  obtained.  [II]  and  Jll]  can  be  easily  illustrated  with  regressions  on  Bocharov’s  published 
yields,  as  explained  earlier  (Figure  19).  There  are  only  a  handful  Balapan  events  with  pub¬ 
lished  yields  in  Bocharov  et  al.  (1989).  However,  the  5  large  historical  events  (for  which  the 
yields  were  exchanged  during  JVE)  can  also  provide  some  further  clue  in  support  of  our  pos¬ 
tulated  hypotheses  [I]  through  [III].  The  yield  estimate  based  on  P^ax  the  three  historical 
events  in  SW  subregion  (790804B,  791223B,  and  811018B)  is  systematically  larger  than  that 
based  on  Pg.  On  the  other  hand,  the  two  events  in  NE  subregion  (791028B  and  840526B) 
have  smaller  yield  estimate  based  on  P^ax  as  compared  to  Pg.  The  larger  bias  of  0.15  m.u. 
that  Ringdal  etal.  (1992)  obtained  with  mt(ISC)  could  have  been  “enhanced”  by  [IV]  and  [V]. 
The  mt,  determination  procedure  presented  in  this  study  does  not  correct  for  [IV]  either.  How¬ 
ever,  the  contribution  of  inter-station  correlation  alone  is  believea  to  be  insignificant  if  WWSSN 
is  used. 

In  Figure  20  we  show  the  difference  of  path  effects  between  BSW  and  BNE  at  each 
WWSSN  station,  which  is  a  measure  of  the  relative  bias  between  BSW  and  BNE  along  each 
path.  Positive  symbols  represent  the  stations  where  BSW  events  are  enhanced  relative  to 
BNE  events.  ISC  network  is  dominated  by  western  European  stations  and  hence  the  effect 
due  to  [V]  would  be  more  severe  than  that  on  WWSSN.  Figure  21  shows  the  spatial  pattern  of 
mi,-Lg  residuals  of  Semipalatinsk  explosions  based  on  Geotech’s  values  and  RMS  Lg 
values  reported  at  NORSAR.  There  is  a  significant  difference  in  the  source  medium  across 
the  Chinrau  fault  separating  the  northeastern  and  southwestern  portion  of  Balapan  test  site, 
as  reported  by  Ringal  and  Marshall  (1989)  and  Marshall  etal.  (1984)  as  noted  in  Table  10. 
The  mean  mi,-Lg  bias  between  SW  and  NE  Balapan  is  about  0.11  m.u.  Figure  21  also  indi¬ 
cates  that  SW  events  near  the  edge  of  the  test  site  tend  to  have  larger  Lg  excitation  (and 
hence  negative  rrib-Lg  residual).  Although  this  seems  to  be  reasonable,  we  must  be  cautious 
as  this  interpretation  is  highly  dependent  on  the  accuracy  of  the  location  as  well  as  the  geo¬ 
logical  information. 
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Top:  RMS  Lg(SAIC)  -  0.782  log(W)  -  4.416 

Bottom:  m2.9(WWSSN)  -  0.768  log(W)  -  4.399 

Lg(top):  RMS  Lg  based  on  Soviet  waveforms  (Israelson,  1991) 

m2.9(bottom):  VWSSN-network  average  of  tmb{l,j)-S{l)-F(l,k(l))] 

Yield  (W):  Bocharov  et  al.  (1989) 

Locations:  Bocharov  et  al.  (1989),  Marshall  et  al.  (1984) 
Tectonics:  Bonham  et  al.  (1980),  Leith  (1987) 


Figure  19.  Regressions  using  yields  published  by  Bocharov  et  al.  (1989)  indicate  that  BNE  explosions 
have  positive  Lg  residuals  (top)  and  negative  mt,  residuals  (bottom);  whereas  BSW  explosions  show  the 
opposite  trend.  Thuc  it  would  seem  plausible  that  the  apparent  mt,-Lg  bias  could  have  been  “enhanced" 
by  the  negative  correlation  between  ivt,  and  Lg  residuals.  There  is  a  distinct  difference  in  the  source 
media  between  the  NE  and  SW  portions  of  Balapan  test  site,  with  the  granites  closer  to  the  surface  and 
the  alluvium  thinner  in  the  southwest.  The  thicker  alluvium  layer  in  NE  region  could  increase  the 
waveform  complexify  and  reduce  the  magnitudes  measured  with  . 
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Figure  20.  Averaged  SW-NE  bias  at  each  WWSSN  station.  Positive  symbols  represent  the  stations 
where  amplitude  of  BSW  events  is  enhanced  relative  to  that  of  BNE  events  of  the  same  source 
strength.  This  pattern  reflects  the  difference  of  path  effects  on  these  two  adjacent  test  sites.  For  net¬ 
work  with  an  uneven  geographical  distribution  of  stations  (such  as  ISC),  the  simple  network  averaging  of 
station  magnitudes  can  only  eliminate  the  path  effect  to  certain  extent. 
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Lg:  NORSAR  RMS  Lg  (RIngdal,  1990;  RIngdal  and  Marshall,  1989) 
Locations:  Bocharov  et  a).  (1989),  Marshall  et  al.  (1984) 

Tectonics:  Bonham  et  al.  (1980),  Leith  (1987) 


Figure  21.  The  spatial  pattern  of  residuals  of  Semipalatinsk  explosions  with  TG's  mi,(GLM)  and 
RMS  Lg  values  reported  at  NORSAR.  The  residual  pattern  of  Balapan  events  strongly  indicates 
significant  difference  in  the  source  medium  across  the  Chinrau  fault  separating  the  northeastern  and 
southwestern  portion  of  the  test  site,  as  reported  by  Ringal  and  Marshall  (1989)  and  Marshall  et  al. 
(1984).  The  mean  m^j-Lg  bias  between  SW  and  NE  Balapan  is  about  0.1 1  m.u. 


52 


6.  PRELIMINARY  ASSESSMENT  OF  WWSSN’S  REMOTE  MONITORING  CAPABILITY 


The  potential  capability  of  a  seismic  network  for  monitoring  a  low  threshold  test  ban 
treaty  is  a  complex  function  of  the  source  and  path  characteristics,  network  geometry,  and  the 
signal-to-noise  threshold  of  the  network's  stations  (Bratt,  1991).  Ringdal  (1990b)  conducted  a 
study  of  the  NORESS  array  detection  capability  for  Semipalatinsk  explosions.  The  50% 
detection  threshold  at  NORESS  is  estimated  as  /nb=3.7±0.1.  Ringdal  noted  a  large  bias  (of 
0.6  mb)  between  Balapan  and  Degelen/Murzhik  subregions.  Thus  the  actual  NORESS  detec¬ 
tion  threshold  at  50%  level  is  mb=2J  and  3.3  for  Balapan  and  Degelen  subregions.  Ringdal 
(1990a)  found  subarrays  of  NORSAR  have  very  different  mj,  ret. 'dual  patterns  for  Novaya 
Zemlya  explosions,  with  mb  bias  ranges  from  +0.9  (03C01)  to  -0.3  (01B05).  For  each 
WWSSN  station,  we  computed  the  sum  of  its  “averaged  station  effect”  and  the  path  term 
associated  with  each  test  site  in  Eurasian  continent.  A  grading  from  A  through  E  is  then 
assigned  to  each  station  according  to  :he  magnitude  of  this  combined  “station  amplification”: 
A  (excellent):  0.3  and  larger;  B  (good):  from  0.1  to  0.3;  C  (fair):  from  -0.1  to  0.1;  D  (accept¬ 
able):  from  -0.3  to  -0.1;  E  (poor):  -0.3  and  smaller  (Table  12). 


Table  12.  WWSSN’s  Capability  in  Monitoring  Eurasian  Explosions 

Station 

Lon 

Lat 

BSW 

BNE 

BTZ 

DEG 

MZK 

AZG* 

NNZ 

PRC 

AAE 

38.766 

9.029 

E 

E 

E 

E 

E 

D 

C 

B 

AAM 

-83.656 

42.300 

A 

A 

A 

C 

D 

A 

A 

B 

ADE 

138.709 

-34.967 

— 

— 

— 

-  - 

— 

— 

- 

B 

AKU 

-18.107 

65.687 

B 

B 

B 

B 

C 

D 

D 

E 

ALQ 

-106.457 

34.943 

— 

— 

— 

... 

— 

B 

E 

f 

1 

1 

ANP 

121.517 

25.183 

E 

— 

E 

C 

E 

— 

E 

C 

AQU 

13.403 

42.354 

D 

D 

D 

D 

D 

E 

A 

B 

ATL 

-84.338 

33.433 

— 

— 

— 

... 

— 

C 

C 

... 

ATU 

23.717 

37.972 

B 

B 

B 

D 

C 

E 

B 

mm 

BAG 

120.580 

16.411 

D 

C 

D 

D 

D 

B 

B 

BEC 

-64.681 

32.379 

C 

D 

C 

D 

E 

D 

D 

... 

*)  AZG:  Azgif  region.  North  Caspian. 
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Station 

Lon 

Lat 

BHP 

-79.558 

8.961 

BKS 

-122.235 

37.877 

BLA 

-80.421 

37.211 

BOZ 

-111.633 

45.600 

BUL 

28.613 

-20.143 

CAR 

-66.928 

10.507 

CHG 

98.977 

18.790 

CMC 

-115.083 

67.833 

COL 

-147.793 

64.900 

COP 

12.433 

55.683 

COR 

-123.303 

44.586 

CTA 

146.254 

-20.088 

DAG 

-18.770 

76.770 

DAL 

-96.784 

32.846 

DAV 

125.575 

7.088 

DUG 

-112.813 

40.195 

EIL 

34.950 

29.550 

EPT 

-106.506 

31.772 

ESK 

-3.205 

55.317 

FLO 

-90.370 

38.802 

FVM 

-90.426 

37.984 

GDH 

-53.533 

69.250 

GEO 

-77.067 

38.900 

GOL 

-105.371 

39.700 

GSC 

-116.805 

35.302 

GUA 

144.912 

13.538 

HKC 

114.172 

22.304 

HLW 

31.342 

29.858 

HNR 

1 59  947 

-9  432 

1ST 

28.996 

41.046 

Table  12.  WWSSN's  Capability  in  Monitoring  Eurasian  Explosions 


BSW  BNE  BTZ  DEG  MZK  AZG 


D 

B 

B 

B 

C 

C 

B 

B 

E 

B 

B 

B 

B 

C 

C 
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Station 

Lon 

Lat 

JCT 

-99.802 

30.479 

JER 

35.197 

31.772 

KBL 

69.043 

34.541 

KBS 

1 1 .924 

78.918 

KEV 

27.007 

69.755 

KIP 

-158.015 

21.423 

KOD 

77.467 

10.233 

KON 

9.598 

59.649 

KRK 

30.062 

69.724 

KTG 

-21.983 

70.417 

LEM 

107.617 

-6.833 

LON 

-121.810 

46.750 

LOR 

3.851 

47.267 

LPS 

-89.162 

14.292 

LUB 

-101.867 

33.583 

MAL 

-4.411 

36.728 

MAN 

121.077 

14.662 

MAT 

138.207 

36.542 

MDS 

-89.760 

43.372 

MSH 

59.588 

36.311 

MSO 

-113.941 

46.829 

MUN 

116.208 

-31.978 

NAl 

36.804 

-1.274 

NAT 

-35.033 

-5.117 

NDI 

77.217 

28.683 

NHA 

109.212 

12.210 

Table  12.  WWSSN’s  Capability  in  Monitoring  Eurasian  Explosions 


BSW  BNE  BTZ  DEG  MZK  AZG 


55 


Station 

Lon 

Lat 

OXF 

-89.409 

34.512 

PDA 

-25.663 

37.747 

PMG 

147.154 

-9.409 

POO 

73.850 

18.533 

PRE 

28.190 

-25.753 

PTO 

-8.602 

41.139 

QUE 

66.950 

30.188 

RAB 

152.170 

-4.191 

RCD 

-103.208 

44.075 

SCP 

-77.865 

40.795 

SDB 

13.572 

-14.926 

SEO 

126.967 

37.567 

SHA 

-88.143 

30.694 

SHI 

52.520 

29.638 

SHK 

132.678 

34.532 

SHL 

91 .883 

25.567 

SJG 

-66.150 

18.112 

SLR 

28.282 

-25.735 

SNG 

100.620 

7.173 

STU 

9.195 

48.772 

TAB 

46.327 

38.068 

TOL 

-4.049 

39.881 

TRI 

13.764 

45.709 

TRN 

-61.403 

10.649 

TUC 

-110.782 

32.310 

UME 

20.237 

63.815 

UNM 

-99.178 

19.329 

Table  12.  WWSSN's  Capability  in  Monitoring  Eurasian  Explosions 


BSW  BNE  BTZ  DEG  MZK  AZG 


C 

C 

D 

B 

B 

B 

C 

C 

B 

C 

C 

E 

E 

D 

E 

D 

E 

E 

C 

C 

E 

B 

C 

B 

C 

c 

D 

C 

B 

C 

C 

D 

E 

AZG 

NNZ 

B 

B 

C 

B 
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Nuttli  (1986)  used  COP  (Copenhagen.  Denmark).  KEV  (Kevo,  Finland),  KON  (Kongs- 
berg.  southern  Nonway).  NUR  (Nurmijarvi.  Finland),  and  UME  (Umea,  Sweden)  in  measuring 
his  mi,{Lg)  values  for  Balapan  explosions.  It  is  interesting  to  note  that  these  five  Scandinavian 
stations  also  give  good  or  excellent  P-wave  recordings.  Nuttli  (1987)  reported  that  none  of 
COP.  KEV.  or  KON  have  readable  Lg  amplitudes  for  the  Degelen  explosions  he  studied, 
although  these  stations  do  have  readable  Lg  amplitudes  for  Balapan  explosions  of  compar¬ 
able  worldwide  network-averaged  m^y.  Our  result  also  show  that  these  three  stations  have 
larger  amplification  for  P  waves  from  Balapan  explosions  as  compared  to  those  from  Degelen 
explosions. 

Butler  and  Ruff  (1980)  examined  SP  P-wave  amplitudes  of  Novaya  Zemlya  explosions 
recorded  by  WWSSN  stations  in  U.S.  The  lowest  amplitudes  were  found  in  GOL  (Golden, 
Colorado)  and  ALQ  (Albuquerque,  New  Mexico),  with  values  a  factor  of  4  lower  than  the  high 
amplitudes.  Our  study  with  an  enlarged  data  set  recorded  at  a  global  network  shows  a  very 
consistent  result.  The  stations  showing  the  lowest  amplitude  are  SHK  (Shiraki,  southern 
Honshu,  Japan),  KTG  (Kap  Tobin,  eastern  Greenland),  GOL,  SJG  (San  Juan,  Puerto  Rico 
region),  LEM  (Lembang,  Java),  and  ALQ.  On  the  other  hand,  COP  (Copenhagen,  Denmark), 
HLW  (Helwan,  United  Arab  Republic),  MSH  (Mashhad,  Iran-Turkmenistan  border),  1ST  (Istan¬ 
bul,  Turkey),  AQU  (Aquila,  central  Italy),  and  ESK  (Eskdalemuir,  United  Kingdom)  report  the 
highest  amplitude  for  NNZ  explosions.  Note  that  the  station  COP,  which  has  the  largest  com¬ 
bined  station  amplification,  is  also  used  in  Nuttli’s  (1986,  1988)  mt,(Lg)  study  of  Balapan  and 
Novaya  Zemlya  events. 
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7.  DISCUSSION  AND  CONCLUSIONS 


Along  with  an  extensive  data  set  of  woridwide  explosions  recorded  at  a  global  network, 
teleseismic  body-wave  amplitudes  from  92  Semipalatinsk  explosions  are  measured  and 
analyzed  to  isolate  the  propagation  characteristics  and  to  derive  a  better  measure  of  the 
source  size.  This  new  mj,  factoring  procedure  reduces  effectively  the  scatter  in  the  station 
magnitudes  across  the  network.  For  the  extreme  case,  the  variation  reduction  could  reach  a 
factor  of  about  3.  In  principle,  it  can  be  applied  to  other  types  of  network-recorded  magnitudes 
as  well,  such  as  mi,  (the  band-passed  spectral  amplitude,  see  Bache,  1982,  and  Murphy 
eta!.,  1989),  my,Lg).  Mo,  and  Ms- 

We  have  recomputed  the  yield  estimates  of  Semipalatinsk  explosions  based  on  the 
path-corrected  values  derived  in  this  study.  Based  on  the  first  motion  of  the  P  wave 
alone,  the  (central-value)  yield  estimate  of  Soviet  JVE  (09/14/88)  is  113  KT  and  that  of  the 
cratering  event  on  01/15/65  at  Balapan  test  site  is  111  KT,  which  are  in  excellent  agreement 
with  estimates  derived  by  other  means.  Thus  the  first  motion  of  the  initial  short-period  P 
waves  appears  to  be  a  very  favorable  source  measure  for  explosions  fired  in  hard  rock  sites 
underlain  by  the  stable  mantle  (such  as  Semipalatinsk).  The  bias  relative  to  NTS  at  50- 
KT  level  is  inferred  as  0.36  magnitude  unit.  Also  included  in  this  study  are  some  preliminary 
assessment  of  WWSSN's  capability  in  remotely  monitoring  Eurasian  explosions.  A  strong 
correlation  between  the  P-wave  amplitude  and  Lg  detection  at  teleseismic  distance  is 
observed. 

Although  our  results  can  explain  some  of  the  propagation  complexities  in  the  initial 
P-wave  arrivals,  a  follow-up  study  is  needed  to  quantify  further  the  contribution  of  near¬ 
source  scattering  to  the  waveform  complexity  in  the  P  coda  which  is  not  covered  in  this  study. 
Our  previous  modeling  effort  (McLaughlin  and  Jih,  1986,  1987,  1988;  Jih  and  McLaughlin, 
1988)  of  utilizing  the  linear  finite-difference  code  (Jih  etal.,  1988)  focused  on  the  effects  of 
mountainous  topography  and  hypothetical  heterogeneity  in  the  upper  crust  on  teleseismic  and 
regional  phases  with  somewhat  simplied  structures  of  other  test  sites.  We  suggest  that  the 
follow-up  research  be  accompanied  with  some  well-constrained  forward  modeling  study  using 
more  realistic  structures  of  ex-Soviet  test  sites. 
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